Isospin-dependent properties of asymmetric nuclear matter in relativistic mean-field 

models 
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Using various relativistic mean-field models, including the nonlinear ones with meson field self- 
interactions, those with density-dependent meson-nucleon couplings, and the point-coupling models 
without meson fields, we have studied the isospin-dependent bulk and single-particle properties of 
asymmetric nuclear matter. In particular, we have determined the density dependence of nuclear 
symmetry energy from these different relativistic mean-field models and compare the results with the 
constraints recently extracted from analyses of experimental data on isospin diffusion and isotopic 
scaling in intermediate-energy heavy ion collisions as well as from measured isotopic dependence 
of the giant monopole resonances in even-A Sn isotopes. Among the 23 parameter sets in the 
relativistic mean-filed model that are commonly used for nuclear structure studies, only a few are 
found to give symmetry energies that are consistent with the empirical constraints. We have also 
studied the nuclear symmetry potential and the isospin-splitting of the nucleon effective mass in 
isospin asymmetric nuclear matter. We find that both the momentum dependence of the nuclear 
symmetry potential at fixed baryon density and the isospin-splitting of the nucleon effective mass in 
neutron-rich nuclear matter depend not only on the nuclear interactions but also on the definition 
of the nucleon optical potential. 

PACS numbers: 21.65.+f, 21.30.Fe, 24.10.Jv 



I. INTRODUCTION 



Besides the many existing radioactive beam facilities 
and their upgrades, many more are being constructed 
or under planning, including the Cooling Storage Ring 
(CSR) facility at HIRFL in China the Radioactive 
Ion Beam (RIB) Factory at RIKEN in Japan @, the 
FAIR/GSI in Germany j|, SPIRAL2 / GANIL in France 
[H, and the Facility for Rare Isotope Beams (FRIB) in 
the USA H. These new facilities offer the possibility 
to study the properties of nuclear matter or nuclei un- 
der the extreme condition of large isospin asymmetry. 
As a result, the study of the isospin degree of freedom 
in nuclear physics has recently attracted much atten- 
tion. The ultimate goal of such study is to extract in- 
formation on the isospin dependence of in-medium nu- 
clear effective interactions as well as the equation of state 
(EOS) of isospin asymmetric nuclear matter, particularly 
its isospin-dependent term or the density dependence of 
the nuclear symmetry energy. This knowledge, especially 
the latter, is important for understanding not only the 
structure of radioactive nuclei, the reaction dynamics in- 
duced by rare isotopes, and the liquid-gas phase transi- 
tion in asymmetric nuclear matter, but also many critical 
issues in astrophysics @, 0, H & EE EH El El- Unfor- 
tunately, the density dependence of the nuclear symme- 
try energy, especially its behavior at high densities, is 
largely unknown and is regarded as the most uncertain 
among all the properties of isospin asymmetric nuclear 
matter. Although the nuclear symmetry energy at nor- 
mal nuclear matter density po rs 0.16 fm~ 3 is known to 



be around 30 MeV from the empirical liquid-drop mass 
formula 14 , |T^] , its values at other densities are poorly 
known [fi ^7 1 . Various microscopic and phenomenological 
models, such as the relativistic Dirac-Brueckner-Hartree- 
Fock (DBHF) 0, El EE El, HE HH, H and the non- 
relativistic Brueckner-Hartree-Fock (BHF) [H, [H] ap- 
proach, the relativistic mean-field (RMF) model based on 
nucleon- meson interactions [T2], and the non-relativistic 
mean-field model based on Skyrme-like interactions [25], 
HE El HE HE EE [HI, have been used to study the 
isospin-dependent properties of asymmetric nuclear mat- 
ter, such as the nuclear symmetry energy, the nuclear 
symmetry potential, the isospin-splitting of nucleon ef- 
fective mass, etc., but the predicted results vary widely. 
In fact, even the sign of the symmetry energy above 
3po is uncertain [32l |. The theoretical uncertainties are 
mainly due to the lack of knowledge about the isospin 
dependence of in-medium nuclear effective interactions 
and the limitations in the techniques for solving the nu- 
clear many-body problem. As to the incomprcssibility of 
asymmetric nuclear matter, it is essentially undetermined 
[33j, even after about 30 years of studies. For compar- 
ison, the incompressibility of symmetric nuclear matter 
at its saturation density po has been determined to be 
231 ± 5 MeV from the nuclear GMR [H[ and the EOS 
at densities of 2po < P < 5po has also been constrained 
by measurements of collective flows in nucleus-nucleus 
collisions 0. 

As a phenomenological approach, the RMF model has 
achieved great success during the last decade in describ- 
ing many nuclear phenomena [IE HE El HE HE HE EL 
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|42[. For example, it provides a novel saturation mecha- 
nism for the nuclear matter, an explanation of the strong 
spin-orbit interaction in finite nuclei, a natural energy 
dependence of the nucleon optical potential, etc. The 
RMF approach is generally based on effective interaction 
Lagrangians that involve nucleon and meson fields. In 
this approach, a number of parameters are adjusted to 
fit the properties of many nuclei. As such, these models 
usually give excellent descriptions of nuclear properties 
around or below the saturation density. 

Since the original Lagrangian proposed by Walecka 
more than 30 years ago [sE], there have been a lot of 
different treatments, extensions, and applications of the 
RMF model. The three main versions arc the nonlinear 
models [H, [3?i [H, Hi| , models with density-dependent 
meson- nucleon couplings [H, 0, 
coupling models without mesons [48 
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For each version of the RMF model, there are also many 
different parameter sets with their values fitted to the 
binding energies and charge radii of a large number of 
nuclei in the periodic table. Including isovector mesons 
in the effective interaction Lagrangians further allows the 
RMF model to describe successfully the properties of nu- 
clei far from the /^-stability line. With recent develop- 
ments in constraining the isospin-dependent properties 
of asymmetric nuclear matter, especially the density de- 
pendence of the nuclear symmetry energy, it is of great 
interest to see to what extend the results from different 
versions of the RMF model are consistent with these con- 
strains. 

In the present work, based on commonly used 23 dif- 
ferent parameter sets in three different versions of the 
RMF model, we carry out a systematic study of the 
isospin-dependent bulk and single-particle properties of 
asymmetric nuclear matter. In particular, we study 
the density dependence of the nuclear symmetry energy 
from these RMF models and compare the results with 
the constraints recently extracted from analyses of the 
isospin diffusion data from heavy-ion collisions based on 
the isospin and momentum-dependent IBUU04 trans- 
port model with in-medium nucleon-nucleon (NN) cross 
sections [HI, [H, HU, the isoscaling analyses of isotope 
ratios in intermediate energy heavy ion collisions [571 ] . 
and measured isotopic dependence of the giant monopole 
resonances (GMR) in even-A Sn isotopes [58|. Among 
these 23 commonly used interactions in nuclear struc- 
ture studies, only a few are found to give symmetry en- 
ergies that are consistent with the empirically extracted 
one. Furthermore, we study the nuclear symmetry poten- 
tial and the isospin-splitting of the nucleon effective mass 
in isospin asymmetric nuclear matter. Our results indi- 
cate that the nuclear symmetry potential at fixed baryon 
density may increase or decrease with increasing nucleon 
momentum depending on the definition of the nucleon 
optical potential and the interactions used. This depen- 
dence is also seen in the isospin-splitting of the nucleon 
effective mass in neutron-rich nuclear matter. In addi- 
tion, the isospin-splitting of the nucleon scalar density in 



neutron-rich nuclear matter is also studied. 

The paper is organized as follows. In Section [Til we 
discuss some isospin-dependent bulk and single-particle 
properties of asymmetric nuclear matter, such as the 
nuclear symmetry energy, the nuclear symmetry poten- 
tial, and the isospin-splitting of nucleon effective mass 
as well as current experimental and/or empirical con- 
straints on these quantities. The theoretical frameworks 
for the different versions of RMF models, i.e., the non- 
linear RMF models, the models with density-dependent 
nucleon-meson coupling, and the nonlinear and density- 
dependent point-coupling models, are briefly reviewed in 
Section [Ml Results on the isospin-dependent properties 
of asymmetric nuclear matter, i.e., the nuclear symmetry 
energy, the nuclear symmetry potential, and the isospin- 
splitting of nucleon effective mass and the nucleon scalar 
densities in neutron-rich nuclear matter, from different 
versions of RMF models are presented and discussed in 
Section IIVI A summary is then given in Section [V] 
For completeness, the isospin- and momentum-dependent 
MDI interaction, which will be used as a reference in some 
cases for comparison, is briefly described in Appendix [A"l 



II. ISOSPIN-DEPENDENT PROPERTIES OF 
ASYMMETRIC NUCLEAR MATTER 

A. Nuclear symmetry energy 

The EOS of isospin asymmetric nuclear matter, given 
by its binding energy per nucleon, can be generally writ- 
ten as 

E(p, a) = E{p, a = 0) + E sym (p)a 2 + 0(a 4 ), (1) 

where p — p n + p p is the baryon density with p n and 
p p denoting the neutron and proton densities, respec- 
tively; a = (p n — p p )/ (pp + p n ) is the isospin asymmetry; 
E(p, a = 0) is the binding energy per nucleon in sym- 
metric nuclear matter, and 



E sym (p) 



1 d 2 E(p, a) 

2 d^ 2 1 



(2) 



is the nuclear symmetry energy. The absence of odd- 
order terms in a in Eq. ([TJ is due to the exchange sym- 
metry between protons and neutrons in nuclear matter 
when one neglects the Coulomb interaction and assumes 
the charge symmetry of nuclear forces. The higher-order 
terms in a are negligible, e.g., the magnitude of the a 4 
term at po is estimated to be less than 1 MeV [59|, HO, [Ml ■ 
Neglecting the contribution from higher-order terms in 
Eq. ([!]) leads to the well-known empirical parabolic law 
for the EOS of asymmetric nuclear matter, which has 
been verified by all many-body theories to date, at least 
for densities up to moderate values. As a good approxi- 
mation, the density-dependent symmetry energy E sym (p) 
can be extracted from E sym (p) ps E(p, a = 1) — E(p, a = 
0), i.e., the energy change per nucleon when all protons 
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in the symmetric nuclear matter are converted to neu- 
trons while keeping the total nuclear density fixed. In 
this sense, the nuclear symmetry energy gives an estima- 
tion of the binding energy difference between the pure 
neutron matter without protons and the symmetric nu- 
clear matter with equal numbers of protons and neutrons. 
It should be mentioned that the possible presence of the 
higher-order terms in a at supra-normal densities can 
significantly modify the proton fraction in /3-equilibrium 
neutron-star matter and the critical density for the direct 
Urea process which can lead to faster cooling of neutron 
stars [6l[63||. 

Around the nuclear matter saturation density po, the 
nuclear symmetry energy E syia (p) can be expanded to 
second-order in density as 



E sym (p) = -Esym(Po) + 



L ( P - Po 



Po 



18 



ym / P ~ PO 



(3) 

where L and A' sym are the slope and curvature parame- 
ters of the nuclear symmetry energy at po, i.e., 



L 



3/Oo 



sym 



9p 



dE sym (p) . 
dp lp 
2 d 2 E sym (p) 



d 2 P 



\P=Po ■ 



(4) 
(5) 



The L and K sym characterize the density dependence 
of the nuclear symmetry energy around normal nuclear 
matter density, and thus carry important information on 
the properties of nuclear symmetry energy at both high 
and low densities. In particular, the slope parameter L 
has been found to correlate linearly with the neutron-skin 
thickness of heavy nuclei and thus can in principle be de- 
termined from measured thickness of the neutron skin of 
such nuclei [H, |M H, S H, S [ZfJ IzJ . Unfortunately, 
because of the large uncertainties in the experimental 
measurements, this has not yet been possible so far. 

At the nuclear matter saturation density and around 
a = 0, the isobaric incompressibility of asymmetric nu- 
clear matter can also be expressed to second-order in a 
as MM 



K(a) « A + A a 



(6) 



where Kq is the incompressibility of symmetric nuclear 
matter at the nuclear matter saturation density and the 
isospin-dependent part [74| 



A 



asy 



A sym 6L 



(7) 



characterizes the density dependence of the nuclear sym- MeV and A a 



asy 

be extracted experimentally by measuring the GMR 
in neutron-rich nuclei. Earlier attempts based on this 
method have given, however, widely different values. For 
example, a value of A asy = —320 ± 180 MeV with a large 
uncertainty was obtained in Ref. [75| from a systematic 
study of the GMR in the isotopic chains of Sn and Sm. 



In this analysis, the value of Kq was found to be 300 ±25 
MeV, which is somewhat larger than the commonly ac- 
cepted value of 230±10 MeV. In a later study, an even less 
stringent constraint of —566 ± 1350 < A asy < 139 ± 1617 
MeV was extracted from the GMR of finite nuclei, de- 
pending on the mass region of nuclei and the number of 
parameters used in parameterizing the incompressibility 
of finite nuclei [33| . Most recently, a much stringent con- 
straint of A'asy = —550 ± 100 MeV has been obtained in 
Ref. [iH from measurements of the isotopic dependence 
of the GMR in even-A Sn isotopes. 

Besides studies of nuclear structure, heavy-ion reac- 
tions, especially those induced by radioactive beams, also 
provide a useful means to investigate in terrestrial labora- 
tories the isospin-dependent properties of asymmetric nu- 
clear matter, particularly the density dependence of the 
nuclear symmetry energy. Indeed, significant progress 
has recently been made both experimentally and the- 
oretically in extracting the information on the behav- 
iors of nuclear symmetry energy at sub-saturation den- 
sity from the isospin diffusion data in heavy-ion collisions 
from the NSCL/MSU [MIIEII1]- Using the isospin and 
momentum-dependent IBUU04 transport model with in- 
medium NN cross sections, the isospin diffusion data were 
found to be consistent with a density-dependent symme- 
try energy of E sym (p) « 31.6(p/pp) 7 with 7 = 0.69 — 1.05 
at subnormal density [55l . |56| , which has led to the ex- 
traction of a value of L = 88 ± 25 MeV for the slope 
parameter of the nuclear symmetry energy at saturation 
density and a value of A asy = —500 ± 50 MeV for the 
isospin-dependent part of the isobaric incompressibility 
of isospin asymmetric nuclear matter [55|, [H, [7(| ■ This 
has further imposed stringent constraints on both the 
parameters in the isospin-dependent nuclear effective in- 
teractions and the neutron skin thickness of heavy nu- 
clei. Among the 21 sets of Skyrme interactions com- 
monly used in nuclear structure studies, only the 4 sets 
SIV, SV, Gcr, and R CT have been found to give sym- 
metry energies that are consistent with above extracted 
one. Using these Skyrme interactions, the neutron-skin 
thickness of heavy nuclei calculated within the Hartree- 
Fock approach is consistent with available experimental 
data [ZO, [7l[ and also that from a relativistic mean-field 
model based on an accurately calibrated parameter set 
that reproduces the GMR in 90 Zr and 208 Pb as well as 
the isovector giant dipole resonance of 208 Pb (76|. The 
extracted symmetry energy further agrees with the sym- 
metry energy E sym (p) = 31.6(/o/po) ' 69 recently obtained 
from the isoscaling analyses of isotope ratios in interme- 
diate energy heavy ion collisions [57l |. which gives L « 65 
1 -453 MeV. The extracted value of 



metry energy. Information on A asy can in principle K. 



-500 ± 50 MeV from the isospin diffusion data is 



also consistent with the value A a 



-550±100 MeV ob- 



tained from recently measured isotopic dependence of the 
GMR in even-A Sn isotopes [HI. We note that the GMR 
only allows us to extract the value of A asy but not that 
of L. These empirically extracted values for L and A sym 
represent the best and most stringent phenomcnological 
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constraints available so far on the nuclear symmetry en- 
ergy at sub-normal densities. Although the behavior of 
the symmetry energy at high densities is presently largely 
undetermined, much of this information is expected to be 
obtained from future high energy radioactive beam facil- 
ities. 



B. Nuclear symmetry potential 

The nuclear symmetry potential refers to the isovector 
part of the nucleon mean-field potential in isospin asym- 
metric nuclear matter. Besides the nuclear density, the 
symmetry potential of a nucleon in nuclear matter also 
depends on the momentum or energy of the nucleon. The 
nuclear symmetry potential is different from the nuclear 
symmetry energy as the latter involves the integration of 
the isospin-dependent mean-field potential of a nucleon 
over its momentum. The nuclear symmetry potential is 
thus a dynamical quantity while the nuclear symmetry 
energy is a thermodynamic quantity, and both are im- 
portant for understanding many physics questions in nu- 
clear physics and astrophysics. Various microscopic and 
phenomenological models have been used to study the 
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28l . |29| . l30l . |31| , and the predicted results vary 
widely as in the case of the nuclear symmetry energy. 
In particular, whereas most models predict a decreasing 
symmetry potential with increasing nucleon momentum 
albeit at different rates, a few nuclear effective interac- 
tions used in some models give an opposite behavior. 

The nuclear symmetry potential was originally defined 
in non-relativistic models. In particular, the nuclear sym- 
metry potential can be evaluated from 



u sym (p,p) 



U n (p,P) ~ U P {p,P) 
2a 



(8) 



where U n (p,p) and U p (p,p) represent, respectively, the 
neutron and proton single-particle or mean-field poten- 
tials In relativistic models, the nuclear symmetry poten- 
tial can be similarly defined by using the non-relativistic 
reduction of the relativistic single- nucleon potentials. 
The nuclear symmetry potential in relativistic models 
therefore depends on the definition of the real part of the 
non-relativistic optical potential or the nucleon mean- 
field potential deduced from the relativistic effective in- 
teractions, which are characterized by Lorentz covari- 
ant nucleon self-energies. In the relativistic mean-field 
approximation, these self-energies appear in the singlc- 
nucleon Dirac equation 



[ 7/i (^-:E£)-(M T + £f)]Vv = 0, 



n,p 



(9) 



as the isospin-dependent nucleon vector self-energy £^ 
and scalar self-energy £f. In the Hartree approximation 
at the static limit, there are no currents in a nucleus or 
nuclear matter, and the spatial vector components van- 
ish and only the time-like component of the vector self- 
energy £° remains. Furthermore, the nucleon self-energy 



is an energy-independent real and local quantity in the 
standard RMF model. 

There are different methods to derive the real part 
of the non-relativistic optical potential based on the 
Dirac equation with Lorentz covariant nucleon vector and 
scalar self-energies. The most popular one is the so-called 
"Schrodinger-cquivalent potential" (SEP). From the nu- 
cleon scalar self-energy T,f and the time-like component 
of the vector self-energy £°, the "Schrodinger-cquivalent 
potential " is given by [77j: 



(10) 

where E'kin is the kinetic energy of a nucleon, i.e., -Ekin = 
E T — M T with E T being its total energy. Eq. pC)]) shows 
that [/sep,t increases linearly with the nucleon energy 
E T or kinetic energy E^ if the nucleon self-energies are 
independent of energy. We note that by construction 
solving the Schrodinger equation with above SEP gives 
same bound-state energy eigenvalues and elastic phase 
shifts as the solution of the upper component of the Dirac 
spinor in the Dirac equation with same nucleon scalar 
self-energy and time-like component of the vector self- 
energy [23| . The above SEP thus best represents the real 
part of the nucleon optical potential in non-relativistic 
models [20l [78j . The corresponding nuclear symmetry 
potential is given by 



U. 



SEP _ J7sEP,ra — UsEP,f 



2n 



(11) 



with a being the isospin asymmetry. 

Another popular alternative for deriving the non- 
relativistic nucleon optical potential in relativistic mod- 
els is to take it as the difference between the total energy 
E T of a nucleon with momentum p in the nuclear medium 
and its energy at the same momentum in free space [7!| , 
i.e., 



C/ OP T,r = E T - v/p 2 + M* 



= Er-yj{E T - £°) 2 - £?(2M T + ££)(12) 

In obtaining the last step in above equation, the disper- 
sion relation 



E T = Y? T + x /p 2 + {Mr + £?) 2 



(13) 



has been used. This definition for the nucleon optical 
potential has also been extensively used in microscopic 
DBHF calculations [§3] and transport models for heavy- 
ion collisions [8l|. For energy-independent nucleon self- 
energies, Uopt, t approaches a constant value of £° when 
\p\ — > oo, unlike the linear increase of L^sep.t with the nu- 
cleon energy. For \p\ = 0, we have Uqpt,t = £f + £? 
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while U S EP.r = + + + S°) 2 /(2M T ). Therefore, 
f/oPT,T displays a more reasonable high energy behavior 
than Usep,t- We note that unlike Usep,t, Uopt,t does 
not give the same bound-state energy eigenvalues and 
elastic phase shifts as the solution of the upper compo- 
nent of the Dirac equation. As in the case of Usep,t, the 
symmetry potential in this approach is defined by 



OPT _ UqPT.ti — Uqpt iP 



sym 



2a 



(14) 



In Ref. [82| . another optical potential was introduced 
based on the second-order Dirac (SOD) equation, and it 
corresponds to multiplying Eq. (fTQ|) by the factor M T /E T , 



at saturation density. As we will see in the following, al- 
though the predicted energy dependence of nuclear sym- 
metry potential at low energy from the RMF models does 
not agree with the empirical Lane potential, it is consis- 
tent with results from microscopic DBHF [l7|, the ex- 
tended BHF with 3-body forces [2J] , and chiral perturba- 
tion theory calculations [88| , which give a Lane potential 
that also stays as a constant or increases slightly with 
momentum for nuclcons with momenta less than about 
250 - 300 MeV/c or with kinetic energies i? k in < but 
decreases with momentum when the momentum is larger 
than about 250 - 300 MeV/c. 



SOD,t 



(E°) 2 ]. (15) 



For energy-independent nuclcon self-energies, Usod,t has 
the same asymptotical value of E° as ?7opt,t when \p\ — > 
oo. For \p\ = 0, we have U SO d,t = E° + ^ + A ^ +Mt E? + 
2{^+s»+M T ) K T 'f) 2 ~ (E?) 2 ]- Thc symmetry potential 
based on the optical potential of Eq. (jX5j) is given by 



u: 



SOD 
sym 



L^SOD.ti — C^SOD,; 

2a 



(16) 



Above discussions thus show that the optical potentials 
defined in Eqs. (Tj"2"|) and (TT5")) have similar high energy 
behaviors, but they may be very different from that de- 
fined in Eq. (flQ|). If we assume that + £° < M T and 
| E^ | ss | S° | , which has been shown to be generally valid 
in the RMF model even at higher baryon densities, we 
have, however, U S ep.t ~ U S od,t ~ U pt,t — Sf + s< ? 
at low momenta (\p\ ps 0), indicating that above three 
definitions for thc optical potential in the RMF model 
behave similarly at low energies. However, it should be 
stressed that, among the three optical potentials defined 
above, only Usep,t is obtained from a well-defined the- 
oretical procedure and is Schrodinger-equivalent while 
C^OPT,r and Usod.t are used here for heuristic reasons as 
they are of practical interest in microscopic DBHF cal- 
culations, transport models for heavy-ion collisions, and 
Dirac phenomenology. 

Empirically, a systematic analysis of a large number of 
nucleon-nucleus scattering experiments and (p,n) charge- 
exchange reactions at beam energies up to about 100 
MeV has shown that the data can be very well described 
by the parametrization U R y m = a — bEy m with a ~ 22 — 34 
MeV and 6 w 0.1 - 0.2 [H, [H, E|. Although the 
uncertainties in both parameters a and b are large, the 
nuclear symmetry potential at nuclear matter saturation 
density, i.e., the Lane potential L^ane [87||, clearly de- 
creases approximately linearly with increasing beam en- 
ergy E'kin- This provides a stringent constraint on the 
low energy behavior of thc nuclear symmetry potential 



Recently, the high energy behavior of the nuclear 
symmetry potential has been studied in the relativis- 
tic impulse (t-p) approximation based on the empiri- 
cal NN scattering amplitude [89( |. The results indicate 
that the nuclear symmetry potential derived from the 
Schrodinger-equivalent potential at a fixed density be- 
comes almost constant when the nucleon kinetic energy 
is greater than about 500 MeV, independent of the pa- 
rameters used in thc analysis. It is further shown that for 
such high energy nucleons the nuclear symmetry poten- 
tial is slightly negative at baryon densities below about 
p = 0.22 fm -3 and then increases almost linearly to pos- 
itive values at high densities. These results provide im- 
portant constraints on the high energy behavior of the 
nuclear symmetry potential in asymmetric nuclear mat- 
ter. Furthermore, with the Love-Franey NN scattering 
amplitude developed by Murdock and Horowitz [9(| HH , 
the intermediate-energy (100 < Eki n < 400 MeV) behav- 
ior of thc nuclear symmetry potential constructed from 
thc Schrodinger-equivalent potential in isospin asymmet- 
ric nuclear matter has also been investigated recently 
[92| . It shows that thc nuclear symmetry potential at 
fixed baryon density decreases with increasing nucleon 
energy. In particular, the nuclear symmetry potential at 
saturation density changes from positive to negative val- 
ues at nuclcon kinetic energy of about 200 MeV. Such 
an energy and density dependence of the nuclear symme- 
try potential is consistent with those from the isospin- 
and momentum-dependent MDI interaction with x = 
(see Appendix [Al for details on this interaction). These 
results thus provide an important consistency check for 
the energy/momentum dependence of the nuclear sym- 
metry potential in asymmetric nuclear matter, particu- 
larly the MDI interaction with x = 0. On the other 
other, the low energy behavior of the nuclear symmetry 
potential at densities away from normal nuclear density is 
presently not known empirically. Experimental determi- 
nation of both the density and momentum dependence of 
the nuclear symmetry potential is thus of great interest, 
and heavy-ion reactions with radioactive beams provides 
a unique tool to extract this information in terrestrial 
laboratories. 
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C. Nucleon effective mass 

Many different definitions for the nucleon effective 
mass can be found in the literature [H [zl]. In the 
present work, we mainly focus on the following three ef- 
fective masses: the Dirac mass M^j (also denoted as 
M* in the present work), the Landau mass Af£ andau , and 
the Lorentz mass M£ orentz . The Dirac mass M^ h&c is de- 
fined through the nucleon scalar self-energy in the Dirac 
equation, i.e., 



Mr 



(17) 



It is directly related to the spin-orbit potential in finite 
nuclei and is thus a genuine relativistic quantity without 
non-relativistic correspondence. We note that the differ- 
ence between the nucleon vector and scalar self-energies 
determines the spin-orbit potential, whereas their sum 
defines the effective single-nucleon potential and is con- 
strained by the nuclear matter binding energy at satu- 
ration density. From the energy spacings between spin- 
orbit partner states in finite nuclei, the constraint 0.55 M 
— -^Dirac — 0-6 M has been obtained on the value of the 
Dirac mass [II III]. 

The Landau mass Af£ andau is defined as M* andaU T = 

Pj§- in terms of the single-particle density of state 
dE T / dp at energy E T and thus characterizes the momen- 
tum dependence of the single-particle potential. In the 
relativistic model, it is given by 93 1 



rir° 

Mwau,r = (£r-£?)(l-^)-(M T 



Since dp/dE T is in principle measurable, the Landau 
mass from the relativistic model should have a compa- 
rable value as that in the non-relativistic model. Em- 
pirically, based on non-relativistic effective interactions 
such as the Skyrmc-type interactions, calculations of the 
ground-state properties and the excitation energies of 
quadrupole giant resonances have shown that a realis- 
tic choice for the nucleon Landau mass is M- 



jM 



Landau 

0.8 ±0.1 [9J,|95J,|96J,|97J. The smaller Landau mass than 
that of nucleon free mass would lead to a smaller level 
density at the Fermi energy and much spreaded single- 
particle levels in finite nuclei [93[ . 

The Lorentz mass Af£ orentz characterizes the energy de- 
pendence of the Schrodingcr-equivalent Potential C/sep,t 
in the relativistic model and is defined as 78] 



A/, 



Lorentz, r 



M T (1 



dUsEP.r ■ 
dE T ■ 



<*-i?>(i-f§> 



dV s 

(M T + H s T )--^+M r -E T 



Mi 



Landau 



dE T 
T + M T - E T 



(19) 

It has been argued in Ref. [78| that it is the Lorentz mass 
-^Lorentz that should be compared with the usual non- 
relativistic nucleon effective mass extracted from analyses 



carried out in the framework of non-relativistic optical 
and shell models. It can be easily seen that in the non- 
relativistic approximation (E T w M T ), the Lorentz mass 
lorentz reduces to the Landau mass M£ andau . 

In relativistic models, the nucleon effective mass has 
sometimes also been introduced via the energy depen- 
dence of the optical potential in Eq. (fT2|) and the second- 
order Dirac optical potential in Eq. (fT5")) . i.e., 



M 



OPT,r 



Mr(1 _ dUo pL) 



AI- 



dE T 

(J5 T -E°)(l-g) 



(^r-Sf)g 



V(£r-£?) 2 -£?(2Af r + £?) 



= M T 



M, 



Landa 



y/(E T - - E?(2A/ T + £#) 



(20) 



and 



M 



SOD,r 



MAl _d_u^ L) 



MJ 



dE T 

Landau, r 



E T 



(Mr + £*) 2 - Cgr ~ Kf + Et-M> 

2EI 



(21) 



respectively. 

The isospin-splitting of nucleon effective mass in asym- 
metric nuclear matter, i.e., the difference between the 
neutron and proton effective masses is currently not 
known empirically [98]. Previous theoretical investiga- 
tions have indicated that most RMF calculations with 
the isovector S meson predict Afjt> iracn < M£ irac while 
in the microscopic DBHF approach, Af£, irac n can be 
larger or smaller than Mp irac _ depending on the ap- 
proximation schemes and methods used for determining 
the Lorentz and isovector structure of the nucleon self- 
energy [20| . For the nucleon Lorentz mass, the micro- 
scopic DBHF or BHF approach and most non-relativistic 
Skyrmc-Hartree-Fock calculations predict M£ orentz n > 
while most RMF and a few Skyrme-Hartree- 



M, 



Lorentz, p > 



Fock calculations give opposite predictions. 



III. RELATIVISTIC MEAN-FIELD MODELS 

For completeness, we briefly introduce in the follow- 
ing the main ingredients in the nonlinear RMF model, 
the density-dependent RMF model, the nonlinear point- 
coupling RMF model, and the density-dependent point- 
coupling RMF model. We neglect the electromagnetic 
field in the following since in the present work we are 
interested in the properties of the infinite nuclear mat- 
ter. Furthermore, besides the mean-field approximation 
in which operators of meson fields arc replaced by their 
expectation values (the fields are thus treated as classi- 
cal c- numbers), we also use the non-sea approximation 
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which neglects the effect due to negative energy states in 
the Dirac sea. 



The nonlinear RMF model 



1. Lagrangian density 

The Lagrangian density in the nonlinear RMF model 
generally includes the nucleon field ip, the isoscalar-scalar 
meson field a, the isoscalar-vector meson field u>, the 
isovector-vector meson field p, and the isovector-scalar 
meson field 5, i.e., 

£nl = $ M 1 ^ - g^) - (M - g a a)] t/> 
+\{d^a - mla 2 ) - \u^ v + \ml^ 

1 



-g P Pfj, ■ ip^fip + gsS ■ i>T^ , 



(22) 



where the antisymmetric field tensors lu^ and p^ are 
given by w M „ = d u uj p - d^ v and p^ v = d u p p - d^p„ , re- 
spectively, and other symbols have their usual meanings. 
Also, vectors in isospin space are denoted by arrows. This 
model also contains cross interactions between the isovec- 
tor meson p and isoscalar a and ui mesons through the 
cross-coupling constants A5 and A v [66|,[9!|. In addition, 
we include the isovector-scalar channel (6 meson) which is 
important for the saturation of asymmetric nuclear mat- 
ter and has also been shown to be an important degree of 
freedom in d escribing the properties of asymmetric nu- 
clear matter |l00l . Il0l| . The above Lagrangian density 
is quite general and allows us to use most of presently 
popular parameter sets in the nonlinear RMF model. 



2. Equation of motion and nucleon self-energies 

From the standard Eulcr-Lagrangc formalism, we can 
deduce from the Lagrangian density equations of motion 
for the nucleon and meson fields. The resulting Dirac 
equation for the nucleon field is 

[ 7|t (i^-S^)-(M + 53f)]^ = > (23) 
with the following nucleon scalar and vector self-energies: 



S t = -9<rO- - gsS ■ f, 
K = g u (j" 



(24) 
(25) 



For the isoscalar meson fields a and lo, they are de- 
scribed by the Klein-Gordon and Proca equations, re- 



spectively, i.e., 
(d^ + m 2 a )a 



g a [4>ip - b a M{g a a) 2 - c a {g„a) 
+A s (g<,<T)g 2 p p li -p»], (26) 



dpuf + mlu" = gSn^-^gl^^u") 
-Avg^pf, ■ p* 1 g u u v ) . 



(27) 



Analogous equations for the isovector d and p meson 
fields are 



(dud* + m 2 5 )6 



gsfrV, (28) 
gpl^Tip - A s {g p p l/ ){g<j(jf 
-k v {g P j?)gluj^]. (29) 



For a static, homogenous infinite nuclear matter, all 
derivative terms drop out and the expectation values of 
space-like components of vector fields vanish (only zero 
components po and luq survive) due to translational in- 
variance and rotational symmetry of the nuclear mat- 
ter. In addition, only the third component of isovector 
fields (5^ and p^) needs to be taken into consideration 
due to the rotational invariance around the third axis in 
the isospin space. In the mean-field approximation, me- 
son fields are replaced by their expectation values, i.e., 
a ^ a, lu^ ^ luq, 5 ^ <S (3) , and p p 
meson field equations are reduced to 



Pg 3 \ and the 



2 -(3) 
m pP0 



9a[ps - b a M(g a a) 2 - c a (g a a) 
+A s {g a a){g p p { £ ) )\ 
gJ[PB - c w (3 w (i> ) 3 

-A(g u u> ){g p p£>) 2 }, 
9s(ps,p ~ PS,n)- 

9p[PB,p ~ PB.n ~ A S (g p p y Q >)(gaVY 

-A v {g p pf ) )(g u u} ) 2 ]. 



(30) 

(31) 
(32) 

(33) 



In the above, the nucleon scalar density ps is defined as 
Ps = = PS.p + Ps,n , (34) 

with the proton (p) and neutron (n) scalar densities given 
by 



PS,, 



(2^r 

Ml 
2tt 2 



dPk 



M* 



k pE p 



(M 



A: 2 

^ 2 ln" 



M* 



where 



E\ 



(M* 



p,n 
(35) 

(36) 



with M* and M* denoting the proton and neutron Dirac 
masses, respectively, i.e., 



M* 
p 



M 



9<yO - gs 



5 {3 \ M*=M- 



9<jO 



gsS^. (37) 
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The nucleon scalar and vector self-energies are then given 
by 



^0 - i -(3) 



(38) 
(39) 



with T3 = 1 and —1 for protons and neutrons, respec- 
tive ly. 



3. Nuclear matter equation of state 

The set of coupled equations for the nucleon and meson 
fields can be solved self-consistently using the iteration 
method, and the properties of the nuclear matter can 
then be obtained from these fields. From the resulting 
energy-momentum tensor, we can calculate the energy 
density e and pressure P of asymmetric nuclear matter, 
and the results are given by 



kin 
1 

+ 2 



kin 
,2 -2 



i 2 -2 , 2r(3)2 , 



2^(3)2 



1 4 3 

+ -b a M{g a a) + -c a [[g a a) + -c u (g u ui ) 

+ l(g p p^f[A s (g aC r) 2 + 3A V { 9u) Q f] (40) 



and 



P 



P, 



kin 

i 

2 



pp 

kin 



2 -2 



mlojl + mjS {3)2 - m-po 



2^(3)2 



--b a M(g a af - -jC a {g a of + -c w {g w u) Q f 
+\{g P p ( i ) ) 2 [Ks{g^f + Avto) 2 ]. 



(41) 



In the above, e l kin and Pl in are, respectively, the kinetic 
contributions to the energy densities and pressure of pro- 
tons and neutrons in nuclear matter, and they are given 
by 



- kin 



(2tt)3 



d^Jk 2 + {M*f 



and 



P 



kin 



-[Z& FPB , + M*ps, 



p,n, 



(42) 



3(2tt) 3 



k 2 



\[E i F p B ,i-M; Ps>i ] 



k 2 + (M*) 2 
i = p, n. 



(43) 



The binding energy per nucleon can be obtained from 
the energy density via 



E 



M, 



PB 



while the symmetry energy is given by 



x 2 x , _ x 2 



E sym (p B ) = i + if^., 



QE F 2 \m 



2 \ms 



M* 2 PB 



E 2 [l+{M-\ A(k F ,M*)} 



(44) 



with the effective p-meson mass given by [66[ 

m f = m l + 9 2 P [As {g a 5-) 2 + Av(g^of} (45) 



and 



A(k F ,M*) 



I"" fc2 
d s k 



W° ffc 2 + (M*) 2 ^ /2 



ps _ Ab 
M* E F 



(46) 



where = \Jk F + M* 2 and M* is the nucleon Dirac 
mass in symmetric nuclear matter. 



B. The density-dependent RMF model 

1. Lagrangian density 

In the density-dependent RMF model, instead of intro- 
ducing terms involving self-interactions of the scalar me- 
son field and cross-interactions of meson fields as in the 
nonlinear RMF model, the coupling constants are den- 
sity dependent. The Lagrangian density in this model is 
generally written as 

-{M -T a a-T s 5-T)]il) 
+ l -(d a ad»a - m 2 s a 2 ) + l -(dJ ■ d»5 - m 2 S 2 ) 
1 .„. 1 



\pau-^ + \m 2 p Pa-(? 1 



(47) 



The symbols used in above equation have their usual 
meanings as in Eq. ([22|) but the coupling constants T a , 
T u , T$ and T p now depend on the (baryon) density, which 
are usually parametrized as 



Flip) = Ti(p sat )hi(x), x = p/ Psat , 



(48) 



with 



l + hjx + d,) 2 . 
hi{x)=ai— — j— — rj, i = a,u>,d,p, (49) 

and p sa t being the saturation density of symmetric nu- 
clear matter. In some parameter sets, 



h p (x) = cxp[-a p (:c - 1)] 
is used for the p meson. 



(50) 
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2. Equation of motion and nucleon self-energies 



3. Nuclear matter equation of state 



Since the coupling constants in the density-dependent 
RMF model depend on the baryon fields and ip through 
the density additional terms besides the usual ones in 
the nonlinear RMF model appear in the field equations 
of motion when the partial derivatives of £dd are car- 
ried out with respect to the fields ip and ip m the Eulcr- 
Lagrange equations. The resulting Dirac equation for the 
nucleon field now reads: 

[ 7| ,(i^-5^)-(M + Ef)]V = 0, (51) 

with the following nucleon scalar and vector self-energies: 

-T a a - T S S '■ f, 



(52) 

(53) 



The new term S^^) in the vector self-energy, which is 
called the rearrangement self-energy (43l l44j , is given by 



P dp 

or 

dp 



dp 



dp 

0T0(5) 



ipTj u i(; ■ p v 



(54) 



with j M = "07^-0. The rearrangement self-energy plays 
an essential role in the applications of the theory since 
it guarantees both the thermodynamical consistency and 
the energy- momentum conservation [43l [44| . 

For the meson fields, the equations of motion are 



ml)a 



(d^+m 2 5 )S = r 5 0T0, 
d v p^ v +m 2 p p^ = Tp^fyj). 



(55) 
(56) 
(57) 
(58) 



In the static case for an infinite nuclear matter, the 
meson equations of motion become 



From the energy-momentum tensor, the energy density 
and pressure of nuclear matter can be derived, and they 
are given by 



t kin ' t kin 
1 



,2 -2 



2 -(3)2 

m P Po 



(66) 



and 



P 



pn 

kin 

1 

~2 



PL 



PB 



2-2 2—2 , 2 T(3)2 

m, a G — + m s d y ' 



2 -(3)2 

l P Po 



.(67) 



It is seen that the rearrangement self-energy does not af- 
fect the energy density but contributes explicitly to the 
pressure. Furthermore, the symmetry energy can be writ- 
ten as 



E S ym(PB) 



hp 
6E F 



PB 



m s 



M* 2 p B 



E F [l+{^) A(k F ,M*)\ 



(68) 



with notations similarly defined as in the nonlinear RMF 
model. 



r p(Pp - Pn), 

m- 5 o^> = T s (ps, P - Ps,n), 
so the nucleon scalar and vector self-energies are 



2 -(3) 

m P Po 



s2 



T 3 , 



-T a a-T s S^ 

r^ + r p ^ 3) r 3 + s w 



with 



^o(B) dT^ _ dV p _( 3 ) 

-W UoPB + -dp- p ° ^- 

dT a _ dT s I(3) , 
—dp- 093 - -d^ 5 {ps * 



Pn) 

~ PS,n)- 



(59) 
(60) 

(61) 

(62) 

(63) 
(64) 



(65) 



C. The nonlinear point-coupling RMF model 



Lagrangian density 



The point-coupling model is defined by a Lagrangian 
density that consists of only nucleon fields. In the present 
study, we use the Lagrangian density of the nonlinear 
point-coupling model of Refs. [48l. [49j . i.e., 



C 



NLPC 



"free 



c 



-If 



••hot 



"dcr 



(69) 
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with 



respectively, with 



^frcc 
£ 4f 



••hot 



£d0 



^(i7„^ - M)V>, (70) 

-\ ctTv{$T~(pil>) ■ (^T7 M V), (71) 
-g/3s(^) 3 - j7sW) 4 
-ITvK^m^^V')] 2 



7tv[W , 'HvV') ' ( - 0T7 M V')]' 



--M^Vh^X^VW) 



<5tv(9 I /^t7 m '0) • (d^^T^ip)- 



(72) 



(73) 



is 



as(#) + Ps(i>i>) 2 + 7sW>V0 3 
-^□W), (77) 
aT S (i>Tip) 

-6 TS 0(ipri>), (78) 

-6 v n(^ip), (79) 

axv(V'^7 A 'V') + 7tv( ? /''^7 m V ; ) ' (^^/^X^'H^VO 
-^□(fn^X (80) 

In the above, □ = d 2 /(c 2 dt 2 — A) denotes the four- 
dimensional d'Alembertian. In the translationally invari- 
ant infinite nuclear matter, all terms involving deriva- 
tive couplings drop out and the spatial components of 
the four-currents also vanish. In terms of the baryon 
density ps and scalar density ps as well as the isospin 
baryon density p 3 = p p — p n and the isospin scalar den- 
sity ps3 = ps,p ~ Ps,n, the nucleon scalar and vector self- 
energies in asymmetric nuclear matter can be rewritten 
as 

= a s p s + /9s Ps + 7s/4 + "TSPS3T3, (81) 
S° = ayp B + 7v/o| + a TV p 3 T 3 + 7TV/O3T3. (82) 



In the above, £ freo is the kinetic term of nucleons and 
C describes the four-fermion interactions while £ hot and 
C dcr contain, respectively, higher-order terms involving 
more than four fermions and derivatives in the nucleon 
field. For the twelve coupling constants in the Lagrangian 
density, as, ay, "ts, "tv, Ps, 7s, 7v, 7tv, <5s, <V, 5 T s, 
and <5tv, the subscripts denote the tensor structure of a 
coupling with "S" "V" and "T" stand for scalar, vector, 
and isovector, respectively. The symbols a;, Si, Pi, and 71 
refer, respectively, to four-fermion or second-order terms, 
derivative couplings, third- and fourth order terms [48], 

Si. 



3. Nuclear matter equation of state 

The energy density e and the pressure P derived from 
the energy-momentum tensor in the nonlinear point- 
coupling RMF model are given by 



' kin 



+ € 



1 



kin 



1 



:"TSPS3 



1 2 1 2 
--a v p + -a TV p 3 

1 



TlSPs 



:7vP 



4 

r7TvP 3 , 



(83) 



2. Equation of motion and nucleon self-energies 



From the variation of the Lagrangian density Eq. 
with respect to ip, we obtain the following Dirac equation 
for the nucleon field: 



[7 M (^-S")-(M + S S )]V = 0, 



(74) 



P = E F p p + E F p n -e p kin -e 



-^ a SPs 



kin 



2 a TSP s3 + ^a v p 



2 a TVP 3 



+ ^Pspt + ^spt + ^7vp 4 + jItvpI (84) 

where E V F and E F are defined as in Eq. |36|) with the 
nucleon Dirac masses 

M* = a s p s + PsPs + IsPs + a T sPS3, (85) 
M* = a s p s + PsPs + ISPs - a T sPS3- (86) 



where the nucleon scalar (£ ) and vector (£ M ) self- 



energies are 



= Vs + Vts-t, 



V + Y!! 



(75) 
(76) 



Furthermore, the symmetry energy in this model can be 
expressed as 



E S ym{PB) 



k F 
6E F 
1 



1 



-OtTYPB 



M* Z PB 



— Otq — 

2 Tb E F [l-a TS A(k F ,M*)} 



(87) 



11 



with notations again similarly defined as in the nonlinear 
RMF model. 



D. The density-dependent point-coupling RMF 
model 



1. Lagrangian density 

For the density-dependent point-couplin g R MF model, 
we use the Lagrangian density of Refs.[52j, |53j, i.e., 



and 



with 



^DDPC — -£frcc + ^4f + ^de 



iii^d» - M)ij), 

1 G s mw)w>) 



(89) 



~2 G T s(p)(tpTip) ■ (ipTljj) 



D s (p){dM){d v U). 



(90) 
(91) 



In the above, £ froe is the kinetic term of the nucleons 
and £ 4f is a four-fermion interaction while £ dcr repre- 
sents derivatives in the nucleon scalar densities. Unlike 
in the nonlinear point-coupling RMF model, the density- 
dependent point-coupling RMF model used here in- 
cludes only second-order interaction terms with density- 
dependent couplings Gi(p) and Di(pi) that are deter- 
mined from finite-density QCD sum rules and in-medium 
chiral perturbation theory [5^, HH • 



2. Equation of motion and nucleon self-energies 



Variation of the Lagrangian Eq. (|88|) with respect to tp 
leads to the singlc-nucleon Dirac equation 

[7^(i^-S' i )-(M + S s )]V = 0, (92) 

with the nucleon scalar and vector self-energies given, 
respectively, by 



= Vs + Vrs-f+^rS 



where 



V s = G s {i>i>) - D 5 n(^), 
V T s = Gts^tiJj), 

V$ = G TV {i>T^i>), 
f) D 

s,5 - — ^(VM^CWO) 



(93) 
(94) 

(95) 
(96) 
(97) 
(98) 

(99) 



fdG s 



= — , 

2 V dp 



dG 



TS 



8Gy 

dp 
dD s 
' dp 



(^V)(^7mV') 



dp 

dGxv 
dp 



{d v {^)){d v {^)) 



(100) 



In the above, we have pu^ = tpj^ip^ where the four- 
velocity u M is defined as (1 — v 2 ) _1 / 2 (l, v) with v being 
the three-velocity vector, and S r s and E£? represent the 
rearrangement contributions resulting from the variation 
of the vertex functionals with respect to the nucleon fields 
in the density operator p. The latter coincides with the 
baryon density in the nuclear matter rest frame. 

In the translationally invariant infinite asymmetric nu- 
clear matter, the nucleon scalar and vector self-energies 
become 



E* 

e° 



GsPs + G T sPS3T~3 (101) 
GvPb + GtvP 3 t 3 + ^ r \ (102) 

with the rearrangement contribution to the self-energy 
<9G 7 



s o (R)= l [ ^ | + ^ |3 



dp 



dp 



dG v 2 



dG T v 2i 
-Pi\ 



dp 



(103) 



Nuclear matter equation of state 



For asymmetric nuclear matter, the energy density e 
and the pressure P derived from the energy-momentum 
tensor in the density-dependent point-coupling RMF 
model are 



-kin 
1 



e lin ~ o G sPs - t; g tsp 2 S3 



-Gyp 1 + -GtvpI, 



(104) 



and 



P = E p F p p + E F p n - e\ 



-G v p 2 + -GtvpI + t; g sPs + ^G T sPs3 



15G S 2 , 

~2-W P3P 
ldG TV 2 



1 dG v 3 

2 dp 9 

. 1 dG TS 2 
f 2^T^' 



where E F and E F arc defined as in Eq. 
effective nucleon masses 

M* = M + G s ps + G T spS3, 
M* = AI + G s ps-G T spS3- 



(105) 
with the 

(106) 
(107) 



As in the density-dependent RMF model, rearrangement 
contributions appear explicitly in the expression for the 
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pressure. Finally, the symmetry energy can be written 
as 

k 2 1 
E sym (p B ) = — ~- + -xGtvPb 

1 M* 2 Pb 
+ ^G T s-^ — (108 

with similar notations as in the nonlinear RMF model. 



pling constants are determined in a self-consistent proce- 
dure that solves the model equations for representative 
nuclei simultaneously in a generalized nonlinear least- 
squares adjustment algorithm (49j . The parameter set 
FKVW for the density-dependent point-coupling model 
arc determined by the constraints derived from the finite- 
density QCD sum rules, in- medium chiral perturbation 
theory, and experimental data of a number of finite nuclei 




IV. RESULTS AND DISCUSSIONS 



A. Nuclear symmetry energy 



Using above models, we have studied the isospin- 
dependent properties of asymmetric nuclear matter. In 
the following, we focus on results regarding the nu- 
clear symmetry energy, the nuclear symmetry potential, 
the isospin-splitting of nucleon effective mass, and the 
isospin-dependent nucleon scalar density in asymmetric 
nuclear matter. For different versions of the RMF model 
considered in the present work, we mainly consider pa- 
rameter sets commonly and successfully used in nuclear 
structure stu dies. In pa rticu lar, we sele ct the parame 
ter sets NL1 [Tol, NL2 [l02|, NL3[lO~ 
FSU-Gold m 



NL-SH 

ha EEcE 



104 

TM1 [1051], PK1 [lpj, FSU-Gold [76|, HA [107[, NLp 
Hp, flOOl fo r the nonlinear RMF model; TW99 

W, DD-M E1 Il08l . DD-ME2 [ToljPKDD [i06| . DD 
93|, DD-F [llOj], and DDRH-cor7j47[ for the density- 
dependent RMF model; and PC-F1 gg, PC-F2_[4g|, PC- 
F3 H|, PC-F4 [H, PC-LA and FKVW gjfor the 
point-coupling RMF model. There are totally 23 param- 
eter sets, and most of them can describe reasonably well 
the binding energies and charge radii of a large number 
of nuclei in the periodic table except the parameter set 
HA, for which to our knowledge there are no calculations 
for finite nuclei. 

We note that all selected parameter sets include the 
isovector-vector channel involving either the isovector- 
vector p meson or the isovector-vector interaction vertices 
in the Lagrangian. The HA parameter set further in- 
cludes the isovector-scalar meson field 6 and fits success- 
fully some results obt ained from the more microscopic 
DBHF approach [l0?t | . The parameter sets NLpS and 
DDRH-corr also include the isovector-scalar meson field 
6, while PC-F2, PC-F4, PC-LA, and FKVW include the 
isovector-scalar interaction vertices. The parameter sets 
NLpS as well as NLp are obtained from fitting the em- 
pirical properties of asymmetric nuclear matter jl00l | and 
describe reasonably well the bindi ng e nergies and charge 
radii of a large number of nuclei [l 1 ll ] . For the DDRH- 
corr, its parameters are determined from the density- 
dependent meson-nucleon vertices extracted from the 
self-energies of asymmetric nuclear matter calculated in 
the microscopic DBHF approach with momentum correc- 
tions, and it reproduces satisfactorily the properties of 
finite nuclei and the EOS from the DBHF approach [47j . 
In the parameter sets PC-F1, PC-F2, PC-F3, PC-F4 and 
PC-LA for the nonlinear point-coupling model, their cou- 



Fig. Q] displays the density dependence of the nu- 
clear symmetry energy E syul (p) for the 23 parameter sets 
in the nonlinear, density-dependent, and point-coupling 
RMF models. For comparison, we also show in Fig. 
Q] results from the phenomenological parametrization of 
the momentum-dependent nuclear mean-field potential 
based on the Gogny effective interaction (25[, i.e., the 
MDI interactions with x = — 1 (open squares) and 
(solid squares), where different x values correspond to 
different density dependence of the nuclear symmetry en- 
ergy but keep other properties of the nuclear EOS the 
same (55| (see Appendix [A] for details). From analyz- 
ing the isospin diffusion data from NSCL/MSU using the 
IBUU04 transport model with in-medium NN cross sec- 
tions, it has been found that the MDI interactions with 
x = — 1 and give, respectively, the upper and lower 
bounds for the stiffness of the nuclear symmetry energy 
at densities up to about 1.2po [55l l5(|. 

It is seen from Fig. [T]that the density dependence of 
symmetry energy varies drastically among different in- 
teractions. In the nonlinear RMF model, while the de- 
pendence on density is almost linear for most parameter 
sets, it is much softer for the parameter sets FSU-Gold 
and HA. The softening of the symmetry energy from the 
latter two parameter sets is due to the mixed isoscalar- 
isovector couplings As and Ay [H, [99[ which modifies 
the density dependence of symmetry energy as seen in 
Eq. (|44|) . For the parameter set NLpJ, it gives a symme- 
try energy that depends linearly on density at low densi- 
ties but becomes stiffer at high densities due to inclusion 
of the isovector-scalar 5 meson. The approximate lin- 
ear density-dependent behavior of the symmetry energy 
for other parameter sets in the nonlinear RMF model 
can also be understood from Eq. (|4"4"| . which shows that 
the symmetry energy at high densities is dominated by 
the potential energy that is proportional to the baryon 
density if the mixed isoscalar-isovector coupling and the 
isovector-scalar d meson are not included in the model. 

The density dependence of the symmetry energy in 
the density-dependent RMF model is essentially deter- 
mined by the density dependence of the coupling con- 
stants T p and T$ of isovector mesons. Most parameter 
sets in this case give similar symmetry energies except 
the parameter sets PKDD and DDRH-corr. Compared 
with other parameter sets in the density-dependent RMF 
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p/p (fm" 3 ) 

FIG. 1: (Color online) Density dependence of the nuclear 
symmetry energy E synl (p) for the parameter sets NL1, NL2, 
NL3, NL-SH, TM1, PK1, FSU-Gold, HA, NLp, and NLp<5 
in the nonlinear RMF model (a); TW99, DD-ME1, DD-ME2, 
PKDD, DD, DD-F, and DDRH-corr in the density-dependent 
RMF model (b); and PC-F1, PC-F2, PC-F3, PC-F4, PC- 
LA, and FKVW in the point-coupling RMF model (c). For 
comparison, results from the MDI interaction with x = — 1 
(open squares) and (solid squares) are also shown. 



model, PKDD gives a very large while DDRH-corr gives 
a very small value for the symmetry energy at saturation 
density. For point-coupling models, all parameter sets 
(PC-F1, PC-F2, PC-F3, PC-F4 and PC-LA) in the non- 
linear point-coupling RMF model predict almost linearly 
density-dependent symmetry energies while the parame- 
ter set FKVW in the density-dependent point-coupling 
RMF model gives a somewhat softer symmetry energy. 

Fig. [TJthus shows that only a few parameter sets can 
give symmetry energies that are consistent with the con- 
straint from the isospin diffusion data in heavy-ion col- 
lisions, which is given by results from the MDI interac- 
tions with x = — 1 and 0. The main reason for this is 
that most parameter sets in the RMF model have satu- 
ration densities and symmetry energies at their satura- 
tion densities which are significantly different from the 
empirical saturation density of 0.16 fm -3 and symmetry 
energy of 31.6 MeV at this saturation density. To show 
this more clearly, we list in Table U the bulk properties 
of nuclear matter at saturation density: the binding en- 
ergy per nucleon —B/A (MeV), the saturation density of 
symmetric nuclear matter po (f m_3 )i the incompressibil- 
ity of symmetric nuclear matter Kq (MeV) , the symmc- 



TABLE I: Bulk properties of nuclear matter at the satura- 
tion point: -B/A (MeV), p (fm" 3 ), K (MeV), E Eym (p ) 
(MeV), K sym (MeV), L (MeV), and K asy (MeV) using the 
23 parameter sets in the nonlinear, density-dependent, and 
point-coupling RMF models. The last column gives the ref- 
erences for corresponding parameter sets. 



Model 


-B/A 


Po 


K 


E Sym 


L 


^sym 


-K a sy 


Ref. 


ATT 1 


16.4 





152 


212 


43.5 


140 


143 


—697 


[102] 


NL2 


17.0 





146 


401 


44.0 


130 


20 


-750 


[102] 


NL3 


16.2 





148 


271 


37.3 


118 


100 


-608 


[103] 


NL-SH 


16.3 





146 


356 


36.1 


114 


80 


-604 


[104] 


TM1 


16.3 





145 


281 


36.8 


111 


34 


-632 


[105] 


PK1 


16.3 





148 


282 


37.6 


116 


55 


-641 


[106] 


FSUGold 


16.3 





148 


229 


32.5 


60 


-52 


-412 


[76] 


HA 


10. D 


n 
U 


170 


233 


30.7 


55 


-135 


-465 


[107] 


NLp 


16.1 





160 


240 


30.3 


85 


■'5 


—507 


[100] 


NLp5 


16.1 





160 


240 


30.7 


103 


127 


-491 


[100] 


TW99 


16.2 





153 


241 


32.8 


55 


-124 


-454 


[46] 


DD-ME1 


16.2 





152 


245 


33.1 


55 


-101 


-431 


[108] 


DD-ME2 


16.1 





152 


251 


32.3 


51 


-87 


-393 


[109] 


PKDD 


16.3 





150 


263 


36.9 


90 


-80 


-620 


[106] 


DD 


16.0 





149 


241 


31.7 


56 


-95 


-431 


[93] 


DD-F 


16.0 





147 


223 


31.6 


56 


-140 


-476 


[110] 


DDRH-corr 


15.6 





180 


281 


26.1 


51 


155 


-151 


[47] 


PC-F1 


16.2 





151 


255 


37.8 


117 


75 


-627 


[49] 


PC-F2 


16.2 





151 


256 


37.6 


116 


65 


-631 


[49] 


PC-F3 


16.2 





151 


256 


38.3 


119 


74 


-640 


[49] 


PC-F4 


16.2 





151 


255 


37.7 


119 


98 


-616 


[49] 


PC-LA 


16.1 





148 


263 


37.2 


108 


-61 


-709 


[49] 


FKVW 


16.2 





149 


379 


33.1 


80 


11 


-469 


[53] 



try energy E sym (p ) (MeV), K sym (MeV), L (MeV), and 
K asy (MeV) using the 23 parameter sets in the nonlin- 
ear, density-dependent, and point-coupling RMF models. 
It is seen that these parameter sets give saturation den- 
sities varying from po = 0.145 fm -3 to po = 0.180 fm -3 
and nuclear symmetry energies -E S ym(po) (MeV) ranging 
from 26.1 to 44.0 MeV. 

To remove the effect due to differences in the satura- 
tion densities among different parameter sets, we show 
in Fig. [5] both the symmetry energy E sym (p) and the 
symmetry energy scaled by its value at corresponding 
saturation density, i.e., E sym (p) / '-E S ym(po) a s functions of 
the scaled baryon density p/po for different parameter 
sets. For comparison, we also plot in the inset in panel 
(b) of Fig. [5] the symmetry energy E sym (p) as a function 
of the baryon density p without scaling. It is seen that 
more parameter sets among the 23 sets become consis- 
tent with the constraint from the isospin diffusion data 
in heavy-ion collisions after scaling the baryon density 
by the saturation density, and with further scaling of the 
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FIG. 2: (Color online) The symmetry energy E aym (p) (a) and 
the scaled symmetry energy E sy m{p) / E synl (po) (b) as func- 
tions of the scaled baryon density p/po for the 23 parameter 
sets in the nonlinear, density-dependent, and point-coupling 
RMF models. Results of the MDI interaction with x = — 1 
(open squares) and (solid squares) are also included for 
comparison. The inset in panel (b) shows the symmetry en- 
ergy E syul (p) as a function of the baryon density p without 
scaling. 

symmetry energy by its value at corresponding satura- 
tion density, most of the parameter sets are in agreement 
with the constraint from the isospin diffusion data. It is 
also interesting to sec from the inset in Fig. [2] that most 
of the parameter sets obtained from fitting the properties 
of finite nuclei give roughly the same value of about 26 
MeV for the nuclear symmetry energy at the same baryon 
density of p = 0.1 fin . This interesting feature is very 
similar to that found with Skyrme interactions [HI, |70| | . 
It implies that the constraint on the symmetry energy 
from fitting the properties of finite nuclei is particularly 
sensitive to the nuclear properties at lower densities, i.e., 
at density slightly above half-saturation density. 

For the density dependence of the nuclear symmetry 
energy around saturation density, a more reasonable and 
physically meaningful comparison is through the values 
of L and A asy given by these parameter sets since the L 
parameter is correlated linearly to the neutron-skin thick- 
ness of finite nuclei while the A asy parameter determines 
the isotopic dependence of the GMR for a fixed clement. 
From Table [H we have seen that the values of L, A sym , 
and lv aS y vary drastically, and they are in the range of 

51 - 140 MeV, -140 - 143 MeV and -750 151 

MeV, respectively. The extracted values of L = 88 ± 25 
MeV and K asy = —500 ± 50 MeV from the isospin dif- 



fusion data, L « 65 MeV and A asy w -453 MeV from 
the isoscaling data, and K asy = —550 ± 100 MeV from 
the isotopic dependence of the GMR in evcn-A Sn iso- 
topes give a rather stringent constraint on the density 
dependence of the nuclear symmetry energy and thus put 
strong constraints on the nuclear effective interactions as 
well. To see this constraint more clearly, we collect in Fig. 
|3]the values of L and A" asy obtained from the 23 param- 
eter sets in the nonlinear, density-dependent, and point- 
coupling RMF models together with the constraints from 
the isospin diffusion data, isoscaling data, and the iso- 
topic dependence of the GMR in even- A Sn isotopes. 
From Fig. [3] as well as Table HI we see clearly that among 
the 23 parameter sets considered here, only six sets, i.e., 
TM1, NLp, NLpS, PKDD, PC-LA, and FKVW, have nu- 
clear symmetry energies that are consistent with the ex- 
tracted L value of 88±25 MeV while fifteen sets, i.e., NL3, 
NL-SH, TM1, PK1, HA, NLp, NLpS, TW99, PKDD, 
DD-F, PC-F1, PC-F2, PC-F3, PC-F4, and FKVW, have 
nuclear symmetry energies that are consistent with the 
extracted A asy value of -500 ± 50 MeV or -550 ± 100 
MeV. Among the latter fifteen sets, only six sets, i.e., 
HA, NLp, NLpS, TW99, DD-F, and FKVW are consis- 
tent with A asy = —500 ± 50 MeV. It is interesting to 
see that most parameter sets in the nonlinear and point- 
coupling RMF models predict stiffer symmetry energies 
(i.e., larger values for the L parameter and larger mag- 
nitudes for A asy ) while those in the density-dependent 
RMF model give softer symmetry energies (i.e., smaller 
values for the L parameter and smaller magnitudes for 

We also see from Table|I]that only five parameter sets, 
i.e., TM1, NLp, NLp<5, PKDD and FKVW, in the 23 
parameter sets have nuclear symmetry energies that are 
consistent with the extracted values for both L and A asy 
(-500 ± 50 MeV or -550 ± 100 MeV). This can be seen 
more clearly in Fig. 0] where the correlation between L 
and A" asy is displayed for the 23 parameter sets together 
with the constraints from the isospin diffusion data, the 
isoscaling data, and the isotopic dependence of the GMR 
in even-A Sn isotopes. Fig. [4] further shows that there 
exists an approximately linear correlation between L and 
A asy , i.e., a larger L leads to a larger magnitude for A asy . 
A similar approximately linear correlation between L and 
A asy has also been observed in Ref. [55| for the phe- 
nomenological MDI interactions, and this correlation can 
be understood from Eq. ([7]) which shows that the value 
of A asy is more sensitive to the value of L than to that 
of A sym . 

The above comparisons thus indicate that the ex- 
tracted values of L = 88 ±25 MeV and A' asy = -500 ±50 
MeV from the isospin diffusion data, L m 65 MeV 
and A asy sa —453 MeV from the isoscaling data, and 
A asy = —550 ± 100 MeV from the isotopic dependence 
of the GMR in even-A Sn isotopes indeed put a very 
stringent constraint on the values of the parameters in 
different RMF models. The fact that most of the 23 
parameter sets considered in the present work give sym- 
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metry energies that are inconsistent with the constraints 
of L = 88 ± 25 MeV and A' asy = -500 ± 50 MeV or 
— 550 ± 100 MeV is probably related to the rather limited 
flexibility in the parametrization of the isovector chan- 
nel in all RMF models. They are also probably con- 
nected to the fact that most of the parameter sets are 
obtained from fitting properties of finite nuclei, which 
are mostly near the /3-stability line and thus are not well 
constrained by the isospin-dependent properties of nu- 
clear EOS. Also, we are interested here in the density- 
dependent behavior of the symmetry energy around sat- 
uration density, as both L and K asy are defined at sat- 
uration density, while the behavior of the nuclear EOS 
at sub-subsaturation density may be more relevant when 
the parameter sets are obtained from fitting the proper- 
ties of finite nuclei. 



B. Nuclear symmetry potential 

Using the parameter sets NL1, NL2, NL3, NL-SH, 
TM1, PK1, FSU-Gold, HA, NLp, and NLp<5 in the non- 
linear RMF model, we have evaluated the energy depen- 
dence of the three different nucleon optical potentials, 
i.e., the "Schrodinger-equivalent potential" U§ep (Eq. 
(fTO)) ). the optical potential from the difference between 
the total energy of a nucleon in nuclear medium and its 
energy at the same momentum in free space Uopt (Eq. 
(fT2")l ). and the optical potential based on the second-order 



Dirac equation Usod (Eq. (|15l )). at a fixed baryon den- 
sity pb = 0.16 fm -3 (roughly corresponding to the satu- 
ration densities obtained from various RMF models). For 
their corresponding symmetry potentials , U®f^ , 

and Ugy^, we have evaluated instead their dependence 
on the nucleon momentum in asymmetric nuclear matter 
at baryon density ps =0.16 fm -3 and with isospin asym- 
metry a = 0.5. We note that in contrast to the energy 
dependence of the nuclear symmetry potential, the mo- 
mentum dependence of the nuclear symmetry potential is 
almost independent of the isospin asymmetry of nuclear 
matter. These results are shown in Fig. [5J Correspond- 
ing results for the parameter sets TW99, DD-ME1, DD- 
ME2, PKDD, DD, DD-F, and DDRH-corr in the density- 
dependent RMF model and for PC-F1, PC-F2, PC-F3, 
PC-F4, PC-LA, and FKVW in the point-coupling RMF 
model are shown in Figs. [5] and [71 respectively. For com- 
parison, we also include in these figures results for the 
energy dependence of the real part of the different opti- 
cal potentials in symmetric nuclear matter at saturation 
density that are extracted from the proton- nucleus scat- 
tering data based on the Dirac phenomenology [H, Ill2l | . 

It is seen that different optical potentials in symmetric 
nuclear matter at pb = 0.16 fm -3 exhibit similar energy 
dependence at low energies but have different behaviors 
at high energies. In particular, at high energies, Usep 
continues to increase linearly with energy while Uopt 
and C/sod seem to saturate at high energies and thus 
display a more satisfactory high-energy limit, similar to 
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what is observed in the nuclear optical potential that is 
extracted from the experimental data based on the Dirac 
phenomenology. The critical energy at which the optical 
potential changes from negative to positive values is be- 
tween about 130 MeV and 270 MeV, depending on the 
parameter sets used. These features are easy to under- 
stand from the fact that the scalar and vector potentials 
are momentum/energy- independent in the RMF models 
considered here. Analysis of experimental data from the 
proton-nucleus scattering in the Dirac phenomenology 
also indicates that the extracted different nucleon op- 
tical potentials in symmetric nuclear matter at normal 
nuclear density change from negative to positive values 
at a nucleon energy of about 208 MeV. Furthermore, it is 
seen that the different optical potentials from all 23 pa- 
rameter sets are consistent with the experimental data at 
lower energies, i.e., Ek in < 100 — 200 MeV, but are gen- 
erally too repulsive at higher energies, especially for the 
"Schrodinger-cquivalent potential" E/sep- These features 
imply that the RMF models with parameters fitted to 
the properties of finite nuclei can only give reasonable de- 
scription of the low energy behavior of the isoscalar opti- 
cal potentials. On the other hand, it should be mentioned 
that for optical potentials at high energies, contributions 
from dispersive processes such as dynamical polarization 
by inelastic excitations, inelastic isobar resonance exci- 
tation above the pi on thresh old, and particle production 
become important [ll3l . Ill4{ . Including such continuum 



> 



CO 

C 

CD 

-I— 1 

O 
Q_ 

TO 
O 

Q. 

o 



200 
100 


-100 
100 


■100 
100 


-100 




Hama et al. 
■ Fit1 

" NLpS(NLp) □ Fit2 





NLpS FSU-Gold NLp - 



NL1 




NL P HA NL P 8 
FSU-Gold 



200 400 600 800 

E-M (MeV) 



j — , — I o 
500 1000 1500 

p (MeV/c) 



FIG. 5: (Color online) Energy dependence of the three differ- 
ent nucleon optical potentials, i.e., Usep (Eq. HQ}), [/opt 
(Eq. flT2j) and [/sod (Eq. (fT5|0 (left panels) as well as 
their corresponding symmetry potentials [/ s ^f , U® 1 ^ , and 
[/ s ^SP as functions of momentum (right panels), at a fixed 
baryon density pB = 0.16 fm -3 for the parameter sets NL1, 
NL2, NL3, NL-SH, TM1, PK1, FSU-Gold, HA, NLp, and 
NLp<5 in the nonlinear RMF model. For comparison, the 
energy dependence of the real part of the optical potential 
in symmetric nuclear matter at saturation density extracted 
from two different fits of the proton-nucleus scattering data 
in the Dirac phenomenology are also included (left panels). 



excitations is expected to improve signi fican tly the high 
energy behavior of the optical potential |l!3l | . Such stud- 
ies are, however, beyond the RMF model based on the 
Hartree level as considered here. 



For the momentum dependence of the symmetry po- 
tential, all 23 parameter sets display similar behaviors in 
U^n, i.e., increasing with momentum, albeit at different 
rates. This can be qualitatively understood as follows. 
Expressing Eq. (fTU)) as 



1 



,?2\ 



(109) 



and neglecting the difference in neutron and proton 
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FIG. 6: (Color online) Same as Fig. [5] for TW99, DD-ME1, 
DD-ME2, PKDD, DD, DD-F, and DDRH-corr in the density- 
dependent RMF models. 



masses, we can rewrite Eq. (fTTj) as 

El - El 



U. 



SEP 



4M T a 
1 



[(E°J 2 + 2E°^ + (A/ n + ES) 2 

S\2 /v0\2 



4M t q 
+ (M„4 

-2E 
1 



° v /p= + (M p + E|)2 - (M p + E*) 2 ] 



[(E^ ) 2 - (E^ 2 + (M^ iraCi „r 



4M T a l 



(110) 



In the simple case of the nonlinear RMF model without 
the isovector-scalar d meson, the neutron Dirac mass is 
the same as that of proton. In this case, U£p£ is reduced 
to 



U. 



SEP 

sym 



0\2 



4M T a 



+2(E°-E«) v /^ + (A/ 5irac ) 2 ]. (Ill) 



Since it can be shown from Eqs. (pH"]) . ( [55]) . and 

that 
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FIG. 7: (Color online) Same as Fig. EJfor PC-F1, PC-F2, PC- 
F3, PC-F4, PC-LA, and FKVW in the point-coupling RMF 
models. 
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(fn - Pp), 



(112) 



we thus have E° > E p and an increase of U s ®^ with the 
momentum of a nucleon in neutron-rich nuclear matter. 
The same argument applies to density-dependent RMF 
models and point-coupling models if the coupling con- 
stant cvtv or Gtv m the point-coupling models is posi- 
tive (at saturation density) so that the potential energy 
part of the symmetry energy at saturation density is also 
positive. 

For f/ S y^ T , whether it increases or deceases with nu- 
cleon momentum depends on the isospin splitting of the 
nucleon scalar self energy (scalar potential) or Dirac mass 
in neutron-rich nuclear matter. This can be seen from 
Eq.(fH|) if it is re-expressed as 



OPT &n Ep 

2^ 



= _(x£-52 + ^ + (M n +Eg)2 



n2 



(113) 

We note that increases with momentum for the 

parameter sets HA, NLp<5, DDRH-corr, and PC-F4 while 
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the opposite behavior is observed for the parameter sets 
PC-F2, PC-LA, and FKVW. 

For the momentum dependence of U^^, it is similar 
to that of if we rewrite Eq. (fTC)) as 



Ul 



SOD 



E„ — E„ 
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E„ 
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E„ 



) 



sym 



4a 



jjOPT 



/2 



(M*+?)(4-- 4-) 



4a 



(114) 



In this case, i/Jy^P increases with nucleon momentum for 
the parameter sets HA, NLpJ, and DDRH-corr while it 
decreases for other parameter sets considered here. 

In Ref . (7^| . it has been argued that it is the 
"Schrodinger-equivalent potential " Usep (Eq. (fTO)) ) and 
thus its corresponding symmetry potential that 
should be compared with the results from non-relativistic 
models. As discussed before, the experimental data in- 
dicate that the nuclear symmetry potential at nuclear 
matter saturation density, i.e., the Lane potential U^ aae , 
clearly decreases at low energies (beam energy E^ n up 
to about 100 MeV and corresponding momentum values 
ranging from about 300 MeV/c to 470 MeV/c), which is 
obviously contradictory to the results for from all 

of the 23 parameter sets considered here. On the other 
hand, and f^fynP f° r some parameter sets can de- 

crease with nucleon momentum, which is qualitatively 
consistent with experimental results. 

For nucleons with momenta less than about 250 — 300 
MeV/c or E^in < 0, although observed increase of 
with momentum for all 23 parameter sets, and f/ S y^ T as 
well as £7 S S >2P with some parameter sets, seems to be con- 
sistent with the results from the microscopic DBHF [l7| , 
the extended BHF with 3-body forces [2J|, and chiral 
perturbation theory calculations [88j . i.e., the symmetry 
potential stays as a constant or slightly increases with 
momentum before decreasing at high momenta, it fails to 
describe the high momentum/energy behaviors of the nu- 
clear symmetry potential extracted from nucleon-nucleus 
scattering experiments and (p,n) charge exchange reac- 
tions at beam energies up to about 100 MeV. 

We note that in studies based on the relativistic im- 
pulse approximation with empirical NN scattering am- 
plitude and the nuclear scalar and vector densities from 
the RMF model, the Schrodinger-equivalent nuclear sym- 
metry potential at fixed baryon density is found to de- 
crease with increasing nucleon energy in the range of 
100 < E ki n < 400 MeV H3 and becomes essentially 
constant once the nucleon kinetic energy is greater than 
about 500 MeV H. 



C. Nucleon effective mass 

For the different nucleon effective masses in symmetric 
nuclear matter at saturation density, we show in Table 
ILTl the results from the 23 parameter sets in the nonlin- 



TABLE II: Values of different nucleon effective masses, 
i.e., M5 irac /M, M L * andau /M, M£ orcntz /A/, M OPT /A/, and 
Msqd/M in symmetric nuclear matter at saturation den- 
sity using the 23 parameter sets in the nonlinear, density- 
dependent, and point-coupling RMF models. The last column 
gives the references for corresponding parameter sets. 



Model 


Mhirac 


M la-ada U 




M OPT 


M son 


Ref. 


M 


Al 


Al 


M 


M 


NL1 




n fid 




n fii 






NL2 


0.67 


0.72 


0.74 


0.70 


0.68 


[102] 


NL3 


0.60 


0.66 


0.67 


0.63 


0.61 


[103] 


NL-SH 


0.60 


0.66 


0.67 


0.63 


0.61 


[104] 


TM1 


0.63 


0.69 


0.71 


0.67 


0.65 


[105] 


PK1 


0.61 


0.66 


0.68 


0.64 


0.62 


[106] 


FSUGold 


0.61 


0.67 


0.69 


0.65 


0.62 


[76] 


HA 


0.68 


0.74 


0.75 


0.71 


0.69 


[107] 


NLp 


0.75 


0.80 


0.82 


0.77 


0.76 


[ 1 001 
lluul 


NLpS 


0.75 


0.80 


0.82 


0.77 


0.76 


[100] 


TW99 


0.55 


0.62 


0.64 


0.60 


0.57 


[46] 


DD-ME1 


0.58 


0.64 


0.66 


0.62 


0.59 


[108] 


DD-ME2 


0.57 


0.63 


0.65 


0.61 


0.59 


[109] 


PKDD 


0.57 


0.63 


0.65 


0.61 


0.59 


[106] 


DD 


0.56 


0.63 


0.64 


0.61 


0.58 


[93] 


DD-F 


0.56 


0.62 


0.64 


0.60 


0.57 


[110] 


DDRH-corr 


0.55 


0.63 


0.64 


0.60 


0.58 


[47] 


PC-F1 


0.61 


0.67 


0.69 


0.64 


0.62 


[49] 


PC-F2 


0.61 


0.67 


0.69 


0.64 


0.62 


[49] 


PC-F3 


0.61 


0.67 


0.69 


0.64 


0.62 


[49] 


PC-F4 


0.61 


0.67 


0.69 


0.64 


0.62 


[49] 


PC-LA 


0.58 


0.64 


0.65 


0.61 


0.59 


[49] 


FKVW 


0.62 


0.68 


0.70 


0.65 


0.63 


[531 



car, density-dependent, and point-coupling RMF mod- 



els. It is seen that the values of M£ irac /M, -^Landau 
Ml olentz /M, M* PT /M, and M* OD /M are in the range 
of 0.55 ~ 0.75, 0.62 ~ 0.80, 0.64 - 0.80, 0.60 - 0.77, 
and 0.57 ~ 0.76, respectively. The parameter sets NL2, 
HA, NLp and NLpi5 seem to give too large values, i.e., 
0.67, 0.68, 0.75, and 0.75, respectively, for the M* irac /M 
as values in the range of 0.55 ~ 0.60 are needed to de- 
scribe reasonably the spin-orbit splitting in finite nuclei 
using the RMF models. On the other hand, the larger 
Dirac masses leads to larger Landau masses M£ andau /M 
of 0.72, 0.74, 0.80, and 0.80, respectively, for the param- 
eter sets NL2, HA, NLp and NLpS, which are consistent 
with the empirical constraint of M* T „j „/M = 0.8 ± 0.1 

The density dependence of the different nucleon effec- 
tive masses in symmetric nuclear matter and correspond- 
ing isospin splitting (M* — M*)/M in asymmetric nu- 
clear matter with isospin asymmetry a = 0.5 are shown 
in Fig. [S] for the parameter sets NL1, NL2, NL3, NL- 
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FIG. 8: (Color online) Density dependence of different 
nucleon effective masses, i.e., M^, i[ac /M , ML andau /M, 
ML orentz /M, Mq PT /M , and Mg OD /M in symmetric nuclear 
matter as well as their corresponding isospin splittings in 
neutron-rich nuclear matter with isospin asymmetry a = 0.5 
for the parameter sets NL1, NL2, NL3, NL-SH, TM1, PK1, 
FSU-Gold, HA, NLp, and NLpS in the nonlinear RMF model. 
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FIG. 9: (Color online) Same as Fig. M but for TW99, DD- 
ME1, DD-ME2, PKDD, DD, DD-F, and DDRH-corr in the 
density-dependent RMF model. 

SH, TM1, PK1, FSU-Gold, HA, NLp, and NhpS in the 
nonlinear RMF model. Figs. 191 and ITTJ1 display the same 
results as in Fig. [5]but for the parameter sets TW99, DD- 
ME1, DD-ME2, PKDD, DD, DD-F, and DDRH-corr in 
the density-dependent RMF models and for PC-F1, PC- 
F2, PC-F3, PC-F4, PC-LA, and FKVW in the point- 
coupling RMF model, respectively. It is seen that dif- 
ferent parameter sets in the nonlinear RMF model give 
significantly different density dependence for the nucleon 
effective masses while the different parameter sets in the 
density-dependent and point-coupling RMF models pre- 
dict roughly the same density dependence for the nucleon 
effective masses except that the parameter set PC-LA 
gives very large values for the nucleon effective masses 
at high densities. This unusual behavior for PC-LA was 
also observed in Ref. (49[, and it is due to the fact that 



FIG. 10: (Color online) Same as Fig. [8] but for PC-F1, PC- 
F2, PC-F3, PC-F4, PC-LA, and FKVW in the point-coupling 
RMF model. 

the coupling constant 7s for the higher-order interaction 
term in PC-LA is positive [48| and dominates at high 
density, leading thus to the very large nucleon effective 
mass. 

For the Landau mass at a fixed baryon density, its 
value Ml andau /M is generally larger than M* irac /Af. 
This can be seen from Eq. (| 18[) if it is rewritten as 

M L * andau>T = (E r - £°) = 0£ _ T + (M T + £?)2 

= y/p 2 F, T + M&Zr (H5) 

which shows that M T * „j „ ^ > MX irar , _ if nucleon self- 
energies are independent of momentum/energy. 

For the Lorentz mass M£ orcntz , M* Lorentz /M depends 
almost linearly on density and thus has a stronger density 
dependence than the Dirac and Landau masses. We note 
from Eqs. (13]) and (fl~T5l) that Eq. (JT9J can be reduced 
to 

Korcntz.r = ~ E?, (116) 

if nucleon self-energies arc independent of momen- 
tum/energy. Therefore, the density dependence of 
M£ orontz is determined uniquely by the density depen- 
dence of nucleon vector self-energy. In the nonlinear 
RMF model, most of the parameter sets, except for TM1, 
PK1 and FSU-Gold which include the self-coupling of 
the to meson field, give a linear density dependence for 
E°, leading thus to the observed linear density depen- 
dence of M£ orontz . As to the nonlinear density depen- 
dence of Af£ orentz in the density-dependent RMF model 
and point-coupling models, it is due to the nonlinear den- 
sity dependence of the coupling constant or the inclusion 
of higher-order couplings. 

For Mq PT /M and M* sou /M, they are seen to have 
roughly same magnitude and also same density depen- 
dence as Af£ andau /M. This feature can be understood 
from the fact that with the dispersion relation of Eq. 
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fT3j) . Eq. (|20|) and Eq. (|21j) can be re-expressed as 



OPT.t 



M£andau.r ( 117 ) 



P F ,t 



and 



SOD,t 



M T 



M 



Landau, r 



(p 2 



A/ 2 ) 



2# 2 



(118) 

respectively. Since p 2 ^ T <C Af 2 (For example, pf ~ 385 
McV/c at p B = 0.5 fm -3 ), we have M T /^p 2 F T + ill 2 « 1 
(with an error of a few percent) and thus A/q PT ss 
M£andauT- Furthermore, the second term in Eq. pi8p 
can be neglected compared with the first term as 
M T /E T ~ 1 (it is a good approximation at low densi- 
ties and with an error of about 20% at high densities, 
e.g., pb = 0.5 fm~ 3 ). As a result, we have MgOD T ~ 
M * 

Landau, r ' 

From the Dirac equation, one sees that condensed 
scalar fields (scalar self- energies) lead to a shift of nu- 
cleon mass such that the nuclear matter is described as a 
system of pseudo-nucleons with masses M* (Dirac mass) 
moving in classical vector fields with 5 meson field or 
isovector-scalar potential further generating the splitting 
of the proton and neutron Dirac masses in asymmet- 
ric nuclear matter. For the isospin splitting of A/ Dirac 
in neutron-rich nuclear matter, it is interesting to see 
that the parameter sets HA, NLp<5, DDRH-corr, and 
PC-F4 give M DiraCjP > M DiraC) „ while PC-F2, PC-LA, 
and FKVW exhibit the opposite behavior of M* Dirac < 
Af£, irac n . This feature implies that the isospin-dependent 
scalar potential can be negative or positive depending on 
the parameter sets used. In the nonlinear RMF model, 
we obtain from Eqs. ([32)) and (|37|) 



M* — M* 

Dirac, n Dirac, p 



-2 ( — ) (PS,n ~ PS,p), 



which indicates that we always have M£> irac > M* Dirac n 
in the neutron-rich nuclear matter where ps, n > 
Ps,p- This argument is also applicable to the density- 
dependent RMF model by replacing g$ with the density 
dependent r<5. For the nonlinear point-coupling models, 
we have, on the other hand, 



M* 

Dirac. n 



-^Dirac,p 



2a Ts(ps,n - PS,p)- 



(120) 



A similar equation can be obtained for the density- 
dependent point-coupling models with the replacement 
of Q!ts by the density dependent Gts- Therefore, the 
isospin splitting of A/ Dirac in neutron-rich nuclear matter 
depends on the sign of the isovector-scalar coupling con- 
stant «ts and Gts in the point-coupling models. Since 
the value of «ts m PC-F2 and PC-LA as well as the value 
of Gts m FKVW are positive, these parameter sets lead 
to the isospin-splitting Mp irac ri > M Dirac in neutron- 
rich nuclear matter, which is opposite to that in other 



parameter sets considered here. The isospin splitting of 
Mpjirac i s directly related to the isovector spin-orbit po- 
tential that determines the isospin-dependent spin-orbit 
splitting in finite nuclei. Unfortunately, there are no clear 
experimental indication about the isospin dependence of 
the spin-orbit splitting in finite nuclei so detailed 
experimental data on the single-particle energy levels in 
exotic nuclei are needed to pin down the isospin splitting 
of M Dirac in asymmetric nuclear matter. 

For the isospin splitting of M£ andau in neutron-rich 
nuclear matter, most parameter sets give M£ andau n > 
M£andau p > which is consistent with the usual re- 
sults in non-relativistic models. The parameter sets 
NLpS and DDRH-corr give, however, the opposite re- 
sult due to the strong isospin-splitting of MQ irac with 
M DiraCi „ < M* Dirac for NLpS and DDRH-corr and the 
fact that M£ andau is related to the Fermi momentum and 
Mfjirac according to Eq. (|115[) . The isospin-splitting 
M £andau,„ > M Lnda.u, P implies that neutrons have a 
larger level density at the Fermi energy and thus more 
compressed single-particle levels in finite nuclei than pro- 
tons. 

For the isospin splitting of M£ orontz in neutron-rich 
nuclear matter, all parameter sets give M£ orentz > 
M£ orcntz ,„ except that the PC-L3 gives M* orentz ' !p < 



M 



Lorcntz,n 



at high densities. From Eq. (|116p . we have 



M 



Lorentz.n 



(121) 



orentz,p 



> 



which leads to the observed isospin-splitting Af£ 
M£ orentz „ as we generally have S° > E° as discussed 
above. For the parameter set PC-L3, it includes a higher- 
order isovector- vector term through the parameter 7tv- 
Since the latter has a negative value and dominates at 
high densities according to Eq. (|52"]) , it leads to S° < S° 
and thus M£ orontz p < M£ orcntz n at high densities. The 
isospin splitting of A/q PT /M and Afg OD /Af in neutron- 
rich nuclear matter show a similar behavior as A/ T * , ,, 
(119) as expected from the discussions below Eqs. (|117p and 



pup . 



D. Nucleon scalar density 



The nucleon scalar density as defined in Eq. is the 
source for the nucleon scalar self-energy (scalar poten- 
tial). In the RMF model, the isospin-dependent nucleon 
scalar density is uniquely related to the nucleon Dirac 
mass as shown in Eq. (|35p . The latter equation also 
shows that the scalar density is less than the baryon den- 



sity due to the factor M* / y k 2 + (A/*) 2 which causes a 
reduction of the contribution of rapidly moving nucleons 
to the scalar source term. This mechanism is responsi- 
ble for nuclear matter saturation in the mean-field theory 
and essentially distinguishes relativistic models from non- 
relativistic ones. In practice, the isospin-dependent nu- 
cleoli scalar density is also an essential ingredient for eval- 
uating the relativistic optical potential for neutrons and 
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protons in the relativistic impulse approximation (See, 
e.g., Refs. [89, 92] and references therein). 
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FIG. 11: (Color online) Neutron and proton scalar densities 
as functions of baryon density in nuclear matter with isospin 
asymmetry a — and 0.5 for the parameter sets NL1, NL2, 
NL3, NL-SH, TM1, PK1, FSU-Gold, HA, NLp, and NLpS 
of the nonlinear RMF model (a); TW99, DD-ME1, DD-ME2, 
PKDD, DD, DD-F, and DDRH-corr of the density-dependent 
RMF model (b); PC-F1, PC-F2, PC-F3, PC-F4, PC-LA, and 
FKVW of the point-coupling RMF model (c). 

In Fig. Ill) we show the neutron and proton scalar 
densities as functions of the baryon density ps in nuclear 
matter with isospin asymmetry a = and 0.5 for the 23 
parameter sets from the nonlinear, density-dependent, 
and point-coupling RMF models. It is seen that the 
neutron scalar density is larger than that of protons in 
neutron-rich nuclear matter at a fixed baryon density. 
Although results for different parameter sets are almost 
the same at lower baryon densities, they become differ- 
ent when ps > 0.25 fm~ 3 , and this is consistent with the 
conclusions of Refs. [H,[92j]. In particular, different pa- 
rameter sets in the nonlinear RMF model predict a larger 
uncertainty for the value of the nucleon scalar density at 
high baryon density while all the parameter sets (except 
PC-LA) in the density-dependent RMF model and point- 
coupling models give roughly same results for the nucleon 
scalar density. These features are consistent with the re- 
sults for the density dependence of nucleon Dirac mass 
shown in Figs. [8j [9j and [TDJ At low baryon densities, 
neutron and proton scalar densities are seen to increase 
roughly linearly with baryon density, and this can be 
easily understood from Eq. (|35p . which is reduced to 



the following expression at low densities (|fc| — > due to 
kp -> 0): 

2 /■*' j3 , M* 

PSi ~ =r / drk — - 

P ' (2tt) 3 Jo M* 

= f F d3k = PB,i> 1 = V, n. (122) 

(27r) Jo 

Therefore, neutron and proton scalar densities generally 
approach their respective baryon densities in asymmetric 
nuclear matter at low baryon densities. 



V. SUMMARY AND CONCLUSIONS 

Using different versions of relativistic mean-field mod- 
els that are commonly used in current nuclear struc- 
ture studies, i.e., the nonlinear model, the model 
with density-dependent nucleon-meson coupling, and the 
point-coupling model, we have investigated systemati- 
cally the isospin-dependent bulk and single-particle prop- 
erties of isospin-asymmetric nuclear matter. In partic- 
ular, we considered 23 parameter sets commonly and 
successfully used in nuclear structure studies, i.e., NL1, 
NL2, NL3, NL-SH, TM1, PK1, FSU-Gold, HA, NLp, 
NLpS for the nonlinear RMF model; TW99, DD-ME1, 
DD-ME2, PKDD, DD, DD-F, and DDRH-corr for the 
density-dependent RMF model; and PC-F1, PC-F2, PC- 
F3, PC-F4, PC-LA, and FKVW for the point-coupling 
RMF model. Most of the parameter sets are obtained 
from fitting the binding energies and charge radii of a 
large number of nuclei in the periodic table or the results 
from the microscopic DBHF approach, which have been 
shown to describe successfully a number of the properties 
of finite nuclei. 

Using these models, we have studied the density de- 
pendence of nuclear symmetry energy and compared the 
results with the symmetry energy recently extracted from 
the analyses of the isospin diffusion data from heavy-ion 
collisions based on an isospin- and momentum-dependent 
transport model with in-medium NN cross sections, the 
isoscaling analyses of isotope ratios in intermediate en- 
ergy heavy ion collisions, and measured isotopic de- 
pendence of the giant monopolc resonances in even-A 
Sn isotopes. These analyses have led to the extrac- 
tion of L = 88 ± 25 MeV for the slope parameter of 
the nuclear symmetry energy at saturation density and 
K asy = -500 ± 50 MeV or -550 ± 100 MeV for the 
isospin-dependent part of the isobaric incomprcssibility 
of isospin asymmetric nuclear matter, which may rep- 
resent the most stringent phenomenological constraints 
available so far on the nuclear symmetry energy at sub- 
saturation densities. Using these constraints, we have 
found that, among the 23 parameter sets considered in 
the present work, only six sets, i.e., TM1, NLp, NLp5, 
PKDD, PC-LA, and FKVW, have nuclear symmetry en- 
ergies that are consistent with the extracted L value of 
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88 ± 25 MeV while fifteen sets, i.e., NL3, NL-SH, TM1, 
PK1, HA, NLp, KLp5, TW99, PKDD, DD-F, PC-F1, 
PC-F2, PC-F3, PC-F4, and FKVW, have nuclear sym- 
metry energies that are consistent with the extracted 
K asy value of -500 ± 50 MeV or -550 ± 100 MeV. Fur- 
thermore, we have found surprisingly that only five pa- 
rameter sets, i.e., TM1, NLp, NLpS, PKDD, and FKVW, 
in the 23 parameter sets have nuclear symmetry energies 
that are consistent with the extracted values for both 
L and K asy . We have noted that most parameter sets 
in the nonlinear and point-coupling RMF models pre- 
dict stiffer symmetry energies while those in the density- 
dependent RMF model give softer symmetry energies. 
These features are probably related to the rather limited 
flexibility in the parametrization of the isovector chan- 
nel in all RMF models and also the fact that most of 
the parameter sets are obtained from fitting properties 
of finite nuclei which are mostly near the /3-stability line 
and thus have less constraint on the isospin-dependent 
properties of asymmetric nuclear matter. Moreover, we 
have focused here on the behavior of the symmetry en- 
ergy around saturation density while the parameter sets 
in RMF models arc fitted to the properties of finite nuclei 
that are more sensitive to the properties of the nuclear 
symmetry energy at sub-saturation densities. 

We have also investigated the energy dependence of 
three different nucleon optical potentials, i.e., the "Schro 
dinger-equivalent potential" C/sep (Eq. (flOj) ). the optical 
potential from the difference between the total energy 
of a nucleon in nuclear medium and its energy at the 
same momentum in free space /7opt (Eq. (|12p ). and the 
optical potential based on the second-order Dirac equa- 
tion J/sod (Eq. ([T5] )), as well as their corresponding 
symmetry potentials , C/ S y^ T , and as func- 

tions of momentum. The results indicate that different 
optical potentials in symmetric nuclear matter exhibit 
similar energy dependence at low energies but have dif- 
ferent high energy behaviors. In particular, at high en- 
ergies, Usep continues to increase linearly with momen- 
tum while C/opt an d C/sod seem to saturate and thus 
display a more satisfactory high-energy limit compared 
to the optical potentials extracted from proton-nucleus 
scatterings using the Dirac phenomenology On the other 
hand, the nuclear symmetry potential at a fixed baryon 
density can increase or decrease with increasing nucleon 
momentum depending on the definition for the nucleon 
optical potential and the interactions used. For 
at pb = 0.16 fm -3 , results from all 23 parameter sets 
show that it increases with momentum, which is consis- 
tent with the predictions of microscopic DBHF and chi- 
ral perturbation calculations at low momenta (less than 
about 300 MeV/c) but is inconsistent with the exper- 
imental result that the nuclear symmetry potential at 
saturation density (the Lane potential l^anc) decreases 
at low energies (beam energy E^in above MeV and less 
than about 100 MeV and corresponding momentum val- 
ues arc from about 300 MeV/c to 470 MeV/c) and RIA 
predictions at higher energies. For and f7~~P , 



they can, however, decrease with momentum for some 
parameter sets, which is qualitatively consistent with the 
experimental constraint. Again, we emphasize that, for 
the three definitions of the optical potential and thus 
their corresponding nuclear symmetry potentials, only 
C^sep,t is well-defined theoretically and is Schrodingcr- 
equivalent while C/opt.t and Usod,t are used here for 
references as £/opt,t has been extensively used in mi- 
croscopic DBHF calculations |8(j and transport models 
for heavy-ion collisions (8l| and U<$ot>,t has been used in 
analyses of the relativistic optical potential based on the 
Dirac phenomenology (82| . 

We have further explored different nucleon effective 
masses, i.e., M } 
M< 



Dirac ' 



g OD m symmetric 



M Landau' M *Lorentz> ^SpT' and 

nuclear matter as well as their 
isospin-splittings in neutron-rich nuclear matter. Most 
of the parameter sets are found to give reasonable values 
for Mp irac as required by the spin-orbit splitting data 
in finite nuclei but too small values for M, 



Landau ' 



imply- 
ing that they would give too small a level density at the 
Fermi energy and too large a spread of the single-particle 
levels in finite nuclei. For M£ orentz , it is found to display 
the strongest (almost linear) density dependence even 
at high densities. Interestingly, including the isovector- 
scalar channel leads to the isospin-splitting of . , and 
-^*Dirac n > -^Dhac p ^ s always obtained in neutron-rich 
nuclear matter for the nonlinear and density-dependent 
RMF models but an opposite result can be observed in 

most parame- 



the point-coupling model. For M-, 



Landau' 

ter sets give the isospin splitting M* an?au „ > M£ andau p 
in neutron-rich nuclear matter, which is consistent with 
usual results in non-relativistic models, while an opposite 
isospin-splitting is observed for M£ orontz . In addition, 
A/q PT , and Mg OU are found to display similar behaviors 



Landa 



as M, 

Finally, we have studied the baryon density depen- 
dence of the nucleon scalar density and its isospin- 
splitting in neutron-rich nuclear matter. The results in- 
dicate that the neutron scalar density is larger than that 
of proton in neutron- rich nuclear matter at a fixed baryon 
density. At low baryon densities, the neutron and pro- 
ton scalar densities generally approach their respective 
baryon densities in asymmetric nuclear matter. 

In the present work, we have focused on three versions 
of standard RMF models, i.e., the nonlinear, density- 
dependent, and point-coupling RMF models. We note 
that there are some recent works [TT5I [TT6I . [TT7I [TT8l [119} 
in which the standard RMF models are extended to in- 
clude density-dependent hadron masses and meson cou- 
pling constants via the Brown- Rho (BR) scaling [l20j |. 
In particular, th e pa r ame ter sets SLC and SLCd con- 
structed in Ref. jll7l . Ill8l ] are not only consistent with 
current experimental results for symmetric matter at nor- 
mal and supra-normal densities and the symmetry en- 
ergy constrained by the isospin diffusion data at sub- 
saturation densities, but also give a fairly satisfactory 
description of the ground state properties of finite nuclei, 
including binding energies, charge radii, and neutron skin 
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thickness. 
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In all standard RMF models, the nucleon self-energies 
arc independent of momentum/energy. As a result, 
the Dirac mass and the Landau mass obtained from 
these models cannot be simultaneously consistent with 
experimental data (see, e.g., Eq. (| 1 1 5[) ) . Also, the 
"Schrodinger-equivalent potential" C/sep (Eq. fTO]) ) in 
these models increases linearly with nucleon energy even 
at high energies. Recently, momentum-dependent nu- 
cleon self-energies have been introduced in the RMF 
model by including in the Lagrangian density the cou- 
plin gs of mes on fields to the derivatives of nucleon densi- 
ties |93l . Il2l| , and the results indicate that a reasonable 
energy dependence of the "Schrodinger-equivalent poten- 
tial" in symmetric nuclear matter at saturation density 
can be obtained, and the Landau mass can also be in- 
creased to a more reasonable value while keeping the 
Dirac mass unchanged, which further leads to an im- 
proved description of /3-decay half-lives of neutron-rich 
nuclei in the Z as 28 and Z m 50 regions (94[. In 
the framework of density-functional theory, including the 
couplings of meson fields to the derivatives of nucleon 
densities in the Lagrangian density provides an effective 
way to take into account higher-order effects. Another 
way to introduce the momentum-dependence in nucleon 
self-energies is to include the Fock exchange terms by 
means of the relativistic Hartree-Fock (RHF) approx- 
imation, even though in practice the inclusion of the 
Fock terms would increase significantly the numerical 
complexity such that it is very difficult to find appro- 
priate effective Lagrangians for the RHF model to give 
satisfactory quantitative description of the nuclear struc- 
ture properties compar e d wi t h standa r d RMF model s 

[M EM EM Ell EM EM EM EH EM, EM EMI. 

Recently, there have been some deve l opm e nts in the 
density-dependent RHF approach jl32l . Il33l . Il34| . It is 
shown that the density-dependent RHF model can de- 
scribe the properties of both finite nuclei and nuclear 
matter with results comparable to those from standard 
RMF models. A more phenomcnological way to improve 
the results of RMF models is to introduce momentum- 
as well as isospin-dependent form factors in the meson- 
nucleon coupli ng constants. It has been shown in Refs. 
[l35l EM Il37| that the empirically observed energy de- 
pendence of the nuclear optical potential in symmet- 
ric nuclear matter at saturation density can be repro- 
duced by relativistic mean-field models with momentum- 
dependent form factors. Finally, to better understand 
the isospin-dependent properties of asymmetric nuclear 
matter it is crucial to investigate the density and momen- 
tum dependence of underlying isovector nuclear effective 
interaction. To reach this ultimate goal, we need not only 
more advanced theoretical approaches but also more ex- 
perimental data both on finite nuclei, especially those 
far from /3-stability line, and from heavy-ion reactions 
induced by high energy neutron-rich nuclei. 
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APPENDIX A: ISOSPIN- AND 
MOMENTUM-DEPENDENT MDI 
INTERACTION 

The isospin- and momentum-dependent MDI interac- 
tion is based on a modified finite-range Gogny effective 
interaction [25[. In the MDI interaction, the potential 
energy density V(p, a) of an asymmetric nuclear matter 
at total density p and isospin asymmetry a is expressed 
as follows [29, ^91 , 



V(p,a) 



Po 



M ( 2 
2/30 



D 



P" 



+ 1 Pi 



X (1 



— Ct.t' 

Po , 



d i pd i p 



3 , fr{r,p)f T ,{r,tf 



l + (p-^) 2 /A 



(Al) 



In the mean-field approximation, Eq. (|Al|) leads to the 
following single-particle potential for a nucleon with mo- 
mentum p* and isospin t in asymmetric nuclear mat- 



U{p,a,p,r) 
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+ B(^) a (l-xa 2 )-8Tx— 

Po a + 



i Pi 



PO 

2C7- — -7- 

Po 



dp 
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i + (p-piy/A 



-ap-i 



(A2) 



In the above r = 1/2 (—1/2) for neutrons (protons); 
(7 = 4/3; f T {r,p) is the phase-space distribution func- 
tion at coordinate f and momentum p. The parameters 
A u (x), Ai(x), B, C T . T , Ct.-t and A are obtained by fitting 
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the momentum-dependence of U(p, a,p,r) to that pre- 
dicted by the Gogny Hartree-Fock and/or the Brueckner- 
Hartree-Fock calculations, the saturation properties of 
symmetric nuclear matter and the symmetry energy of 
31.6 MeV at normal nuclear matter density po = 0.16 
fm -3 [25|. The incompressibility Kq of cold symmetric 
nuclear matter at saturation density po is set to be 211 
MeV. The parameters A u (x) and Ai (x) depend on the x 
parameter according to 

2_£> IB 
AJx) = -95.98 -x , AAx) = -120.57 + x — — . 

cr + 1 cr + 1 

(A3) 

The different x values in the MDI interaction are in- 
troduced to vary the density dependence of the nuclear 
symmetry energy while keeping other properties of the 



nuclear equation of state fixed [55j, and they can be ad- 
justed to mimic the predictions of microscopic and/or 
phenomenological many-body theories on the density 
dependence of nuclear matter symmetry energy. The 
last two terms in Eq. (|A2[) contain the momentum- 
dependence of the single-particle potential. The mo- 
mentum dependence of the symmetry potential stems 
from the different interaction strength parameters C T ._ r 
and C T)T for a nucleon of isospin r interacting, respec- 
tively, with unlike and like nucleons in the background 
fields. More specifically, we use C r ,_ r = —103.4 MeV 
and C T , r = -11.7 MeV. 

With f T (r,p) = tj0(p/(t) — p) for nuclear matter at 
zero temperature, the integrals in Eqs. (|A1|) and (|A2[) 
can be calculated analytically and we find 
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l + (p-2?) 2 /A 2 
A 2 { Pf (r)p f (r') [3(p}(r) +p}(r')) A 2 



+4A 



(pj(r)-pj(r')) tan 



-i P f( T ) -Pfi T ') 
A 



(p1(t)+ P %t')) tan 



-i p/O) +Pf(i j y 

A 



+ i [A 4 + 6A 2 (p 2 (t) +p) { r>)) - 3(p}(r) - p) { r>)f] In + + ^ 



(A4) 



,3 / fr(r,f) 

dp 



2 A3 



P 2 f (r) 



A 2 



2pA 



V ln (p + p/(r)) 2 



A 2 



(P-P/(T))2+A2 



A 



2 tan 



-i P + Pfir) 
A 



2 tan 



A 



(A5) 



With above results as well as the well-known contribution 
from nucleon kinetic energies in free Fermi gas model, 
we can thus easily obtain the EOS of asymmetric nuclear 
matter at zero temperature. 

We note that the MDI interaction has been extensively 
used in the transport model for studying isospin effects 
in intermediate energy heavy-ion collisions induced by 



neutron-rich nuclei [H, El IM EH, [Hi, QH EM EM 

Il43| and the study of th e ther mal properties of asymmet- 
ric nuclear matter jl44l . Il45| . In particular, the isospin 
diffusion data from NSCL/MSU have constrained the 
value of x to between and —1 for nuclear matter den- 
sities less than about 1.2po (55L l56|. 
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